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SUMMARY 


In the fan stage of a turbofan engine, broadband noise is generated by turbulent flow impinging 
on rotor and stator leading edges and by turbulence passing over trailing edges. This report 
provides an analytical model for prediction of the leading edge sources at both the rotor and 
stator. Whereas previous modeling treated rotors and stators as isolated blade rows (thereby 
ignoring reflection/transmission effects of neighboring blade rows), the analysis herein addresses 
broadband noise generation in a fully coupled environment. Modal scattering by the stator and 
mode/frequency scattering by the rotor are included in the acoustic/vortical coupling. This brings 
in the mode trapping physics previously included only in analyses of tone noise. Rotors and 
stators are treated as rectilinear cascades immersed in a mean flow that is constant in the 
spanwise direction. The cascade response theory recognizes 3D perturbations so that a true 3D 
turbulence spectrum can be used for excitation. Flow turning at both blade rows is handled by 
unsteady actuator disk theory based on the 4 linearized equations for conservation of mass and 
momentum. 

If each acoustic mode/frequency combination were coupled to all others, the coupling analysis 
would be unmanageable. However, in developing the analysis, an important principle was 
discovered. The set of all participating modes divides itself naturally into many mode sub-sets 
that only couple within themselves. This property of “independent mode sub-sets” makes the 
broadband problem tractable and is the basis for the analysis. It was found that simple modal 
averaging techniques can be applied to eliminate many of these sub-sets from calculations and 
thereby save considerable computational effort. 

This report presents the derivation of the coupling theory in detail and provides documentation 
for the associated computer code CupBB. Computed results show that accounting for coupling to 
adjacent blade rows adds considerably to downstream noise and increases the predicted split 
(or differential) between upstream and downstream sound power. Effects are large enough to 
indicate that absolute level predictions in the future should include the coupling physics. 

The broadband coupling method developed here can be adapted to CFD procedures in the future. 
Lessons learned regarding modal sub-sets and modal averaging will carry over with little 
modification. 
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SECTION 1 
INTRODUCTION 


Fan noise is generated when flow non-uniformities interact with rotor blades and stator vanes as 
suggested in Figure 1. Sources can be classified into leading or trailing edge types. Leading edge 
sources include rotor inflow disturbances such as inlet distortion, atmospheric turbulence, and duct 
boundary layer turbulence. Stator inflow disturbances are primarily from the rotor wake and endwall 
flow. Trailing edge noise (also called self noise) is generated when turbulence in a blade or vane 
boundary layer interacts with a trailing edge. Fan noise can be further classified into tones or 
broadband as indicated in Figure 2. Tone noise is generated when the rotor interacts with a fixed flow 
distortion and when the stator interacts with the periodic component of the rotor wakes. This report is 
concerned with broadband noise caused by turbulence. The focus herein is on leading edge sources 
although the coupling methodology could be extended to trailing edge sources (self noise) as well. 




Fgure 1. Noise generation in a turbofan by Figure 2. Fan noise spectrum showing major 

turbulence in rotor wakes impinging on stator. contributing sources. 

Modeling for turbulence leading edge noise dates back at least to the early 1960’s when simple 
estimates were made of overall noise based on turbulence intensity. No spectral distributions could be 
made with those early models. In the ensuing years, broadband modeling has become increasingly 
sophisticated with inclusion of spectral shapes, non-compactness, 3 dimensionality, and cascade 
effects. However, all of the preceding modeling treated noise generation by an isolated blade row 
only. For example, in noise generation at the stator in Figure 1, the only role of the rotor would be to 
generate turbulent wakes. Effects of acoustic reflection and transmission were ignored in earlier work. 


The purpose of this report is to provide a first look at the effects a neighboring blade row on broadband 
noise generation. For noise generation at the stator, rotor reflection and transmission loss are 
obviously issues to consider. However, the report goes beyond that by including full unsteady 
coupling of the pressure and vortical waves between a rotor and stator. In addition to reflection and 
transmission losses at both blade rows, this accounts for mode and frequency scattering and mode 
trapping. 


This report is an extension of previous work by the author: Reference 1 treats coupled rotor/stator tone 
interaction in 2D; References 2 and 3 analyze broadband noise generated by isolated rotors and stators. 
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The present report combines these 2 methodologies in a coupled broadband model that has been coded 
as CupBB. This work has been enabled by the appearance of a harmonic cascade theory by Glegg 
(ref. 4) that this report adapts to the broadband problem. Glegg’s theory accounts for inflow 
disturbances with wavenumber in all 3 coordinate directions; hence, 3D turbulence can be represented 
rigorously. 

In the following, Section 2 presents background on the tone theory for coupled cascades and on the 
broadband theory for isolated cascades. It then outlines the strategy for combining these into the 
coupled broadband model. Section 3 summarizes the equations to be solved for the mean and unsteady 
flow. Section 4 develops the concept of “independent mode sub-sets” that facilitates the analysis. 
Sections 5 and 6 derive the coupling equations. Section 7 solves the coupled equations for the 
broadband application and shows how they are driven by a turbulence spectrum, Section 8 provides 
some verification of the new theory and presents sample calculations of sound power spectra to 
establish the importance of coupling effects. Section 9 gives some concluding remarks. Appendices A 
and B show how Glegg’s theory is adapted to the current application. Appendix C derives equations 
for unsteady actuator disks used to represent flow turning at the rotor and stator. Finally, Appendix D 
provides a list of notation. 
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SECTION 2 

BACKGROUND & STRATEGY FOR COUPLED BROADBAND THEORY 


This section provides brief reviews of earlier theories for tone noise generation in a coupled 
rotor/stator environment and for broadband noise from an isolated rotor or stator. Then the strategy is 
outlined for combining these methodologies into the coupled broadband prediction scheme. 


Review of Coupled Tone Method 

Reference 1 presents a 2D scheme for coupling a rotor and stator modeled by flat plates as in Figure 3. 
The rotor sends acoustic and vortical waves to the stator and the stator sends acoustic waves to the 
rotor in a full unsteady coupling. Excitation of the system was via the periodic component of the rotor 
wake. It was found for certain vane/blade 
ratios that the blade passing frequency 
component of the wake produced high levels 
of 2xBPF and 3xBPF noise through 
frequency scattering at the rotor and 
trapping of the fundamental interaction 
mode between rotor and stator. After 
developing the coupling scheme via 2D 
modeling, it was extended to quasi 3D in 
Ref. 5 and is used at Pratt & Whitney for fan 
tone prediction in TFaNS, the Theoretical 
Fan Noise Prediction System. We will 
outline operation of the system in a 
2 dimensional context but readers interested 
in details should consult Refs. 1 and 5. The 
coupling scheme is rederived later in this 
report for broadband application so that this 
report does not depend on the references. 

Waves are considered in 3 regions per Figure 3 and counted by the index r. Also, there are 3 wave 
types: 7=1, 2 for upstream/downstream propagating pressure waves and 7 = 3 for vortical waves in 2D. 
The pressure waveform for a type 7 wave in region r is written in a Fourier series as 

rf (*.') = *>„£ I T = 1,2 (2.1) 

n k 

Coordinate x corresponds to the fan axis and 0 is used for the transverse coordinate. Circumferential 
mode order given by the classic Tyler-Sofrin mode selection rule (Ref. 6): 

m = nB\ — (2.2) 


OR 

A 


REGION 1 



REGION 3 


ROTOR 


STATOR 


Figure 3. Flat plate representation for coupled rotor/ 
stator tone noise in Ref. 1 . 
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Bj and B 2 are the numbers of blades and vanes. The first summation is over the blade passing 
harmonics, counted by n. The second summation is over k, the stator scattering index. Vortical modes 
are represented by their transverse velocity component in parallel form 

vf (x,0 = fl o XX Aj(n,k)e ilk < T » k - y+m<l> ~" B ' nt] , T = 3 (2.3) 

n k 

Scattering by the blade rows is denoted by coefficients defined in the following notation that relates 
modal amplitudes to each other 

Ar’(n, k f ) <= Sf r (n\ k'; n, k)A f (n, k ) (2.4) 

/ 

where Sjf(n',k';n,k) is the scattering coefficient giving the ratio of output wave amplitudes (with 
the primes) to input wave amplitudes (without primes). When all of these scattering coefficients are 
found from a cascade acoustic theory, the resulting equations are arranged in matrix form as follows 
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This partially compressed notation suppresses the n, k subscripts shown explicitly in Eq. 2.4. Dots in 
the matrix represent zero blocks. The notation can be further compressed so that Equation 2.5 reads 

A = SA + B (2.6) 

B is the source vector or array of coefficients of prescribed waves that excite the system. Solution of 
this coupled system for the state vector A is given formally by 

A = (l-Sr 1 B (2.7) 

For wake excitation, the only non-zero elements of B would be the B^(n,k) 's to represent vortical 
waves ( T = 3) from region #2. These coefficients, in an expression like Eq. 2.3 would define the rotor 
wake (only the k = 0 term is needed for convected disturbances). 
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When all of the elements of the state vector are found from Eq. 2.7, the entire wave field can be 
constructed from Equations 2.1 and 2.3. 

Over a range of vane/blade ratios, the fundamental interaction mode can be trapped between rotor and 
stator. In this case, the tone coupling effect can be very strong and actually dominate noise generation. 
For example. Fig. 4, reproduced from Ref. 1 , shows the downstream sound power in 3 BPF harmonics 
that results from exciting the stator with the BPF fundamental of the rotor wake in a coupled 
environment. Below BPF cuton in the aft duct, which is at M, = 0.80, the noise is dominated by energy 
scattered up in frequency by the rotor. We now know how to avoid this situation by proper choice of 
vane count (and thereby, mode selection). However, since broadband noise generation includes all 
modes, mode trapping could also be a factor in that spectrum component and may be more difficult to 
control. One objective of the model developed this report is to determine the importance of mode 
trapping in broadband fan noise. 



Rotor Rotational Mach Number 


Figure 4. Coupled tone calculation from Ref. 1 showing 3 BPF harmonics resulting from excitation of the stator 
by the one harmonic of the wake at BPF. 
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Review of Broadband Method for Noise from an Isolated Stator 

Here we review how the harmonic acoustic response of a cascade can be used in a broadband model 
for noise caused by inflow turbulence. This and the above coupled harmonic analysis will be used to 
explain how the coupled broadband formulation will be derived. 


First, consider the cascade harmonic problem 
with mean flow and upwash as in Figure 5. We 
write the harmonic upwash waveform as 

Vl’(X, /) = Wq / ( W-+#+V---a>f) (2.8) 

In this notation, (f> = y/R represents angle around 
an “unwrapped annulus” of effective radius R. 
z takes the role of the radial coordinate and 
v is its associated wavenumber. With axial 
wavenumber 


*v3<7 = 



(2.9) 


y 



Figure 5. Sketch for cascade discussion. 


Eq. 2.8 represents a convected wave (as in the frozen gust model of turbulence). Say we have a 
cascade response function F k {(0,v) from an analysis like Glegg’s (Ref. 4) that gives the upstream 
acoustic pressure resulting from the upwash as 

p(x,t) = WqY, F k {(0,v)J lkx ^ xHq ~ kBl * +v =- m] (2.10) 

k 


where k is the scattering index for the stator and k xlqk is the axial wavenumber for acoustic waves. 
With these equations, the harmonic acoustic pressure field can be computed from the inflow 
perturbation. 

So far the discussion has related to harmonic excitation and response; w q is a Fourier coefficient and 
Eq. 2.10 is a Fourier series. Now, we want to extend the above for application to broadband noise, in 
which case we deal with Fourier transforms and Fourier integrals. We generalize the inflow 
representation of Eq. 2.8 to aperiodic waves 


m'(x,0 = X JJ H.(«>,v)e' ( W- m) d(odv (2.11) 

<7 

which represents the upwash disturbance as a sum/integral over circumferential order q , radial order 
v , and frequency (0 . (In a more general type of flow, an integral over axial wav enumber would also 
be included. However, for convected waves, axial wavenumber and frequency are related through 
Eq. 2.9 and the fourth integral is not required.) There is an issue of existence of the Fourier transform 
Wq in Eq. 2.11. However, there are standard methods to deal with this, as applied later in the formal 
derivations. 
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With inflow represented by Eq. 2.1 1, the corresponding aperiodic response (generalizing Eq. 2.10) is 


P(x,t) = X X jj »’ q ( 0 >y)F k (ay)e ,{t ''* xHq - l ^ ) * w '--‘‘ a dmdv 

k q 


( 2 . 12 ) 


Note that w q {(0, v ) is now a Fourier transform rather than a Fourier coefficient but that F k ( CO, V') is the 
same acoustic harmonic response function used for the tone analysis. 

In principle, any non-periodic excitation and response fields could be represented in these forms. This 
is not actually possible in the case of turbulent flow, however, because m’(x,/) [and vi’^(a>, v)] cannot be 
known in sufficient detail. Instead we deal with statistical quantities (or expected values). In the 
following we review how the expected value of the pressure mean square (and its spectrum) can be 
related to the turbulence spectrum in an isolated cascade analysis. This is provided as a point of 
departure for the coupled broadband analysis. 

To start, we square Eq. 2.12 and take the expected values of both sides, giving 


l^p 2 (x,t)^ = X X XX 11 11{ w q { 61 , v) M *'(«', v')J> F k {( 0 ,V) Fp{(o',v') 


k k q q 

Xe i\(k x i q k -k xlq ' k ')x+[(q-kB 2 )-(q'-k'B 2 M+(v-v'):-{a)-co')t} ^ J v ' 

(2.13) 


where ( ) is the expected value operator. It can be shown that the expected value of the product of 
velocity transforms is related to the upwash component of the turbulence spectrum as follows 


(u,’ q (oi,v) w*'{(o',v')^ = ~ < 5 ^7 < 7 , $ (co- co')8(v -v')® ww (K) 


(2.14) 


Turbulence wavenumber K has components 


K = (k 



(2.15) 


The delta functions in Eq. 2.14 eliminate 2 integrals and a summation, reducing Eq. 2.13 to 


( P 2 <M))=^X S JJ«>„ w (K >/=-*(«,>•) (®. V ) rftu ( 2. 1 6) 


(k-k')B 2 6] 


This result for the mean square pressure varies with position via the <j) variable. If we limit our interest 
to the average over <p, the averaging process produces another delta, 8 kk r , and we arrive at our final 
form for the mean square pressure integrated over all frequencies 

(p 2 (x,t)) = -j^Y< SjJ°M-HXK)|F^(6),v)| 2 i/ft)Jv (2.17) 

k q 
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where it can be seen that the x dependence has disappeared. (Actually, this is true only for the 
propagating waves.) Usually, we want the sound spectrum or power spectral density of the pressure; 
this is just the frequency integrand: 

S pp ((0) = ^JO vn ,,(K)|^.(ft),v)| 2 dv (2.18) 

k q 

This establishes a very simple relation between the spectrum of turbulence and the sound spectrum 
through the square of the magnitude of the acoustic response function. There is no phase information; 
only magnitudes appear. We are generally not interested in phase information for the sound spectrum 
and this is fortunate since turbulence models do not include phase. 

Figure 6 compares computed results from this theory with scaled data from an ADP fan model. In 
applying the 2D geometry and mean flow, conditions at the 85% radius were used. Turbulence 
intensity and scale were the only free parameters and they were adjusted for a good fit to the data. 
Further detail is given in Ref. 2. The fit is reasonably satisfactory: overall shape, including high 
frequency roll off, is good and the upstream/downstream split matches that of the data. However, it 
will be seen that the match to the upstream/downstream split was fortuitous since the more complete 
analysis shows that the rotor blocks much of the upstream noise from the stator. 



Figure 6. Calculations from BBCascade compared with scaled test data. Plot reproduced from Ref. 2. 
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Broadband Coupling 

Now we ask how to combine the 2 types of analysis above into a coupled broadband model. The key 
is Eq. 2.18: think of the turbulence spectrum O w , M ,(K) as excitation for the system and think of 

as the system response function. There is a parallel with the harmonic solution found in 

Eqs. 2.1 and 2.7; in this case the system response function (1-S)' 1 is the inverse of a large matrix of 
Fourier coefficients. In Section 6 we will see that (1-S)' 1 takes the role of F*(et), v) in Eq. 2.10. An 
integral/summation like Eq. 2.18 will be wrapped around this entire matrix system. In the broadband 
analysis, the elements of S will be Fourier transforms instead of Fourier coefficients. This approach 
retains all of the phase information in the frequency and mode scattering between the blade rows but 
then gives up the phase information in the final working formulas, as was the case with Eq. 2.18. 
Again, the sound spectrum will be driven by the turbulence spectrum. 

The analysis is simpler than might be expected because of the property of “independent mode sets” 
mentioned in the Summary and Introduction. Sections 3 to 6 establish the coupling matrix system, 
Section 7 connects the matrix system to the turbulence spectrum and gives the final working equations 
for the computer code CupBB, and Section 8 shows some computed examples. 
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SECTION 3 

MEAN FLOW AND GENERAL WAVE CHARACTERISTICS 


In this section we define the environment for the coupled blade row analysis. This includes the duct 
boundaries, mean flow, and boundary conditions. It also gives the general properties of the 
perturbation waves consistent with this ducted flow environment (but without yet accounting for the 
presence of blades and vanes). 

In later sections, we use Glegg’s unsteady cascade analysis (Ref. 4) for the acoustic and vortical 
responses of the blade rows. Since this treats waves with 3 independent wavenumbers, it permits a 
completely general representation of turbulent flow. However, it is based on rectilinear geometry and 
a uniform background flow. Because of this, we represent flow in a fan duct annulus by “unwrapping” 
it and enforcing periodic boundary conditions in the tangential direction (i.e., variables at <p =±7t are the 
same). We introduce an “effective radius” R, which is used for scaling and to make the equations look 
more like the true 3D flow equations. R relates the tangential angular and linear coordinates through 

y = <PR- 

Geometry and Mean Flow 

The .y, y, : coordinates are shown in Figure 7. ,v is parallel to the machine axis, y is the tangential 
coordinate (parallel to the direction of rotor rotation), and r corresponds to the radial coordinate. Inner 
and outer duct walls are represented by the r = 0 and z = h planes. Periodic 




Figure 7. Definition of boundaries and mean flow 


boundaries are at y = ±nR. Mean flow is the same everywhere with axial component U and tangential 
(or swirl) component V. Represented in vector form, the mean velocity is 

U = £/i+V j + Ok (3.1) 
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Perturbation Equations 

Here we follow the general approach and notation of Smith (Ref. 7) but extend it to 3 dimensions. 
Unsteady flow satisfies the linearized continuity and momentum equations 


where the convective derivative is 


-^ + p,.a;V- u = 0 (3.2) 

P,^- = ~VP (3-3) 

Dt 

ir,=T, +v - v (34 > 


and subscript r implies mean value in a region designated by r. For perturbations periodic in space 
and time, solutions have the general form 


p 


p 
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u 
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V 
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i(,k x x+k y 


y+k-z-ojt) 


Note that the convective derivative yields 

— = i(-co + U k x +Vk v ) 

Dr A y 

and, accordingly, define a parameter that occurs frequently in the analysis 

A = -©+ Uk x +Vk y 


(3.5) 


(3.6) 


(3.7) 


When we substitute the solutions from Eq. 3.5 into the continuity and momentum equations, the 
exponentials cancel and we can write 


A 

a r Pr k . 

a 2 p r k y 

a 2 p r k. 
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k x /p r 
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0 
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ky/pr 
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_k.Jp r 
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w 


Non-trivial solutions are obtained from 

A 2 [A 2 - a 2 (k 2 + k 2 + k ; )] = 0 (3.9) 
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This permits 4 wave types. For periodicity in y, we require m to be an integer in k y y = m<j>. This, with 
0 =y/R, gives k y = m/R. Also, for hard wall boundary conditions, we use linear combinations of 
Eqs. 3.5 to give sines and cosines for r dependence. Thus, ju is an integer and k z = {fin !h). The 
following form fits all 4 wave types 


p T = p T e KK-r^-co,) cos( Mi) 
u T =u T e^kx-fX+mip-cot) cos( Ml) 

h (3.10) 

v T = v T e i(kxr cos(^) T = 1, 2, 3, 4 

w T =w T A^-^sin(Hi) 

These, with the correct relationships between the pressure and velocity magnitudes (given below), 
satisfy the governing equations, the periodicity conditions in the y direction, and zero radial velocity at 
the duct walls at r = 0 and h. Characteristics of the different wave types will be discussed below. 


Pressure Waves 

In Eq. 3.9 pressure waves are associated with 


[X 2 -a 2 (k 2 +k 2 +kt)] = 0 


This can be solved for the axial wavenumber k x : 


k x = 


of. -U 2 L 


U(-0) +Vk y ) ± a.^i-a+Vkyf -{a 2 -U 2 )(k 2 y + kl) 


(3.11) 


(3.12) 


As in Smith’s analysis, we identify upstream-going pressure (or acoustic) waves as Type 1 waves: 

k xX =- T -^r(/(-a) + wV//?)-c7,. A /(-a) + my /)?) 2 -(a, 2 -t/ 2 )[( w >^) 2 +( M %') 2 ] (3.13) 

at — U~ L 


and downstream-going pressure waves as Type 2 waves: 


k x2 =- T ±-^\u(-(0 + mV/R) + a r J(-<0 + mV/R) 2 -(a 2 -U 2 )[{%) 2 +e n / h ) 2 } 
at -U~ L 


(3.14) 


Depending on frequency and mode orders m and /u, the argument of the square root can be positive or 
negative. If it is positive, the axial wavenumber is pure real and waves propagate without decay. If the 
argument is negative, the axial wavenumber is complex and the signs of the square roots must be 
chosen so that waves decay exponentially away from the blade row in the x direction. 
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For the pressure waves we use the momentum equation to relate the complex velocity magnitudes to 
the complex magnitude defining the pressure waves: 


u T = 


v t — 


w T — 


A' vt — 

m/R - 

L t-i 


Pr 'V 


Pt 


- 


Pr 


Pt 


T = 1,2 


(3.15) 


Vortical Waves 

In Eq. 3.9 vortical waves are associated with A = 0. i.e., 

~(Q + Uk x +V k y = 0 (3.16) 

which can be solved for their axial wavenumber. Since A was squared, there are 2 independent vortical 
waves with the same wavenumber 


_ c o-mV/R 
~ U 


T = 3,4 


(3.17) 


Substitution into the momentum equation confirms that the vortical waves have no pressure 

p T = 0 T = 3,4 (3.18) 

Also, substitution of Eq. 3.17 into the exponentials of Eq. 3.10 verifies that the type 3 and 4 waves are 
purely convected. The relationships between the velocity components depend on continuity and the 
actual families to be used in the modeling. These will be defined in Section 5. 
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SECTION 4 

PERMITTED MODES AND THE SCATTERING RULES 


We intend to build a rotor/stator coupling scheme based on matching modal output of the stator to 
modal input of the rotor and vice versa. The usual tool for studying modal behavior in the fan 
environment is the “mode order/frequency plot” schematically shown at the right in Figure 8. In our 
convention, modes with positive m are co-rotating (rotating in the same direction as the rotor). Such 
plots permit us to locate mode orders at 


Frequency 



Figure 8. Mode order/frequency plot for n = 0 


various frequencies and to determine whether 
they propagate (are cut on) or decay (cut off). 

Modes in the shaded area are cut on; hence, 
we can see that that there are more modes 
propagating at higher frequencies. The 
cutoff boundaries are given by the square 
root in Eqs. 3.13 and 3.14. When the 
argument is positive, modes propagate. Test 
data shown in Figure 9 were plotted in the 
same format. Acoustic power can be found 
at all integer mode orders in the cuton area. 

When the radial modes are included, we find 
that thousands of modes must be accounted 
for to get a complete representation of the 
noise field. 

We will see that, when a mode impinges on 
the stator, it scatters out a series of modes of 
different circumferential order. When each 
of these impinges on the rotor, scattering 
there produces modes of different orders and 
frequencies. If any mode could scatter into 
any other mode and frequency, the modal 
accounting scheme would be unmanageable. 

We will soon see that this is not the case but 
that only modes within certain sub-sets of the 
total mode count are connected with each other; there is no scattering from one sub-set to another. 
Once this is understood, the fan broadband analysis can be set up as a series of “small”, independent 
coupling problems and the situation becomes manageable. This section is devoted to explaining the 
kinematics of broadband scattering/coupling and establishing the independent modal sub-sets just 
mentioned. 


10000 


7500 


5000 


2500 



0 - 

-40 *30 *20 -10 0 10 20 30 40 

Spinning Order m (positive co-rotating) 

Figure 9. Experimental mode order-frequency plot from 
Boeing’s 18 inch research fan, Ref. 8. 


Permitted Modes 

Modal scattering should already be understood in the tone context to readers familiar with the Tyler- 
Sofrin fan noise theory (Ref. 6). In the tone case, the frequencies and modes are restricted to 
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frequency : co = nB^O . , n = ... - 1, 0, 1, 2, ... 

mode order : m = nB\ - kB 2 k = 1, 0, 1, 2,... 


i.e., frequencies are restricted to multiples of blade passing frequency (counted by n) and 
circumferential mode orders are restricted to the values of m resulting from integer values of the rotor 
and stator scattering indices, n and k. Furthermore, it was shown in Ref. 1 that, even after scattering 
back and forth between rotor and stator, waves are still restricted the frequencies and mode orders in 
Eq. 4.1. Thus, Eq. 4.1 is a mode set that can scatter only among its members but not to other orders 
and frequencies. Broadband scattering is more general: any frequency and any integer mode order can 
occur. In the following, we will find how Eq. 4.1 is generalized for the broadband problem. 

The space-time dependence of the pressure waves in Equation 3.10 is given by 

e Hk x .v+m<p-cor) cos(j u^^/ } (4.2) 


where we can identify circumferential mode order m and radial mode order fd. A major 
simplification in the analysis arises when we note that there is no scattering from one radial mode to 
another. This is an artifact of using rectilinear geometry and constant mean flow and is true because 
the pressure modes and upwash components of the vortical modes (see Eq. 3.10) have cosine 
dependence on the radial (r) coordinate. A cosine input produces a cosine output with the same fd. 
This simplification would not apply strictly in a situation with real geometry and flow but, even there, 
the radial coupling should be weak and may not be a major factor in noise generation. Thus, radial 
order fd becomes a parameter and, for coupling at any x plane, we can focus on circumferential orders 

via the exponential e i ^ m< k~ a> 0 j 0 establish the scattering rules, we define nuft - cot to be the kinematic 
phase 


y/ -m(j)-0)t 


(4.3) 


and track its behavior as waves scatter back and forth between the rotor and stator. 


To start the discussion, we must define interblade phase angle <7 and show what happens to it upon 
scattering by a cascade. Since a is the phase shift from y = 0 to y = g (the blade gap), it follows 
(using the exponential form in Eq. 3.5) that a = k y g. In our rectilinear geometry the tangential 
wavenumber k y = m/R. Cascade theory (see Ref. 4 or 7) teaches us that, for an input wave with 
frequency ft) and circumferential wave number k y =o/g , the scattered waves have the same frequency 
ft) and their circumferential wavenumbers are 




<7 - Ink 


(4.4) 


where k is an integer that we call the scattering index and the prime indicates a scattered wave. 
Consider a wave input to the stator with circumferential order q and frequency C0q. The 
circumferential wavenumber of the input wave is k y = <j/g = q/R. And since B 2 g-2nR , it follows that 
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(7 = 


(4.5) 


2nq 

"sT 

And for k' y = m/R , it can be seen from Eq. 4.4 that the orders of the scattered modes are given by 

m = q- kB 2 (4.6) 

Thus, from Eq. 4.3 the kinematic phases of the waves scattered by the stator with input y/ = q(p-0) o t 
are 


y/ — (q — kB 2 )<p — (0 o t (4.7) 

To consider scattering at the rotor, shift to rotor coordinates via 

(j) = (p r + kit (4.8) 

Then, the same mode viewed in the rotor frame has phase 

yr' = (q- kB 2 )0 r - [c 0 o -{q- kB 2 )Q]t (4.9) 

It can be seen that the frequency in that frame is 

c o =CQ 0 -(q- kB 2 )£2 (4. 10) 

For scattering by the rotor, we use the index n (with the opposite sign for convenience) so that the 
waves scattered by the rotor have kinematic phase 


i f/"=(q + nBy - kB 2 )<p r - [co 0 -{q- kB 2 )ft]f (4.11) 

Scattering does not change frequency in the rotor frame. When this expression is shifted back to the 
stator frame via 0 r =0 -kit, we find that the same waves have the form 

yt" = (q + nB\ - kB 2 )<j> - (0) o + wBjQ) t (4.12) 

Mode orders have been shifted from the input order q by nB\ - kBz and frequencies have been 
scattered up and down from the input frequency (0 0 by multiples of blade passing frequency nB/Q . 

This process of following the scattering from stator to rotor and back can be repeated indefinitely but it 
turns out that the waves, when represented in the stator frame, always have the form given by Eq. 4.12. 
Furthermore, we could have started with a wave impinging on the rotor and found the same form. 
Thus, to summarize, a wave in the stator frame having frequency tub and mode order q scatters as 
follows 
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frequency : co 0 => 

mode order : q => 


co = (0 o + nB { Q., 


(4.13) 


n = ...- 1 , 0 , 1 , 2 ,... 
m = q + nB i -kB 1 k =...- 1,0, 1,2,... 

This result establishes the scattering “independent mode sub-set” property asserted above: if we input 
a set of modes with frequencies co 0 +nB ] £l and orders m = q + nB l -kB 2 , those modes only scatter 

among themselves and not to other frequencies and not to other mode orders. In particular, at any 
frequency, only one out of every Bj modes is involved in the scattering. 

We can consider co 0 and q to be frequency and mode offsets from the tone interaction set of Eq. 4.1. 
The set of modes and frequencies resulting from a choice of jjl, and q is independent of sets with 
other values of these indices. Thus, in the computer code, we can set up loops on ji, and q and 
solve a relatively small coupled problem for scattering on n and k for each of these combinations. In 
each case, the input is a vector of turbulence modes for all n and k and the response is found 
simultaneously for the acoustic modes for the same ranges of n and k. Fortunately, the required 
range of k is limited by acoustic cutoff. The range of n to be included could be limited by the 
frequency content of the turbulence; however, practically speaking, this range may be limited by 
computer resources. What we hope is that the n (frequency) range can be truncated at a reasonable 
value without significantly sacrificing accuracy for the frequencies included in the analysis. This did 
turn out to be the case for the tone problem (Ref. 1) and will be verified later in the broadband case. 

This independent mode sub-set property (including the absence of coupling between radial orders) 
changes the broadband coupling problem from one that appeared intractable to one quite manageable, 
at least in terms of the rectilinear geometry used herein. The operation of the independent sub-set 
property will be made clearer below with illustrations in terms of mode order/frequency plots. 


Mode Order/Frequency Plots 

The plot shown in Figure 9 was constructed for a fan with 18 blades and 45 vanes and applies to the 
axial flow regions upstream of the rotor and downstream of the stator. In general, the cuton boundaries 
can be determined by setting the argument of the root in Equation 3.13 to zero: 

(-db + mMy) 2 -p 2 (m 2 +V 2 ) = 0 (4.14) 

where d) = (0R/a r and v = finR/h . This can be re-arranged as follows 

fi 2 co 2 = (coM v + p 2 nt) 2 + ft 2 ft 2 v 2 (4.15) 

which is the equation of a conical surface enclosing the cut on values of m and v for any frequency. 
(To visualize the cone, recognize in Figure 9 that the cone axis is the go coordinate, the v axis is out of 
the paper, and that the plot is a cut through the v =0 plane.) 
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Frequency in Stator Frame 

/V 

(0 



Figure 9. Mode order/frequency plot for fan in axial flow at M x = 0.5. 5/ = 18, Z?2 = 45, and Mj = 0.5. Also, q 
and (do = BPF/2. Presented in stator reference frame. 


Frequency in Rotor Frame 
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In the v =0 plane, Eq. 4.16 reduces to equations for straight lines with slopes M y ± fi x . 

d> = (M y ± P x )m , v = 0 (4.16) 

In Figure 9, the slopes are ±/? v since M x = 0. The large dots on the plot represent permitted modes 
according to Eq. 4.13, for a particular values of co 0 and q. Scattering at the stator is along the dashed 
lines; this is at constant frequency (constant n ) and varying k. Scattering at the rotor is along the 
dotted lines; this is at constant k and varying n. It is easily deduced that these lines have slope M T 
and that this slope is less than that of the co-rotating cutoff boundary, Af v +/?, , for cases where the flow 
is subsonic relative to the rotor blades. 

Figure 10 has exactly the same information as Figure 9 but presented in the rotor reference frame. 
Frequency at the rotor is given by Eq. 4.10. Changing reference frames does not change what modes 
are cut on; it only changes their apparent frequencies. 

Behavior of the independent modal sub-set is now established in Figure 9: with the stator scattering 
along horizontal lines and the rotor scattering along the diagonal lines, any mode inside the cuton 
limits can communicate with any other such mode through multiple reflections between rotor and 
stator. However, only those modes and frequencies are connected with each other. Figure 9 was 
plotted with q = 0 and the broadband frequencies half way between the BPF harmonics (0J o = 
BiQH). Other independent mode sets can be found by taking the dashed and dotted scattering lines as 
a rigid pattern and moving the pattern vertically on the page to cover different frequencies and 
horizontally to cover different q' s. Each new position (each new q and (0 o ) generates a new set of 
modes that scatters only among themselves. Furthermore, for radial orders greater than zero, the 
pattern of dashed and dotted scattering lines can be held above the plane of Figure 9 to represent 
cutting the conical surface for jj = 1, 2, 3, ... for more independent mode sets. Note, with this example 
for 45 vanes, q can take on all integer values from 0 to 44. The broadband frequency variable is 
continuous so that any number of tub’s could be used; however, a broadband spectrum could be 
reasonably well defined by frequencies at the BPF multiples and halfway in between. 

The plots just discussed were generated for axial flow, which would apply upstream of the rotor and 
downstream of the stator Between blade rows, swirl affects the cuton boundaries as shown in 
Figure 11. Both the swirl flow and axial flow boundaries are plotted. The gray region in the plot 
indicates modes that propagate between rotor and stator but are cut off upstream and downstream. 
These modes are “trapped”; that is they can bounce back and forth between blade rows, possibly 
amplifying, but cannot escape. However, they do scatter to higher mode orders and frequencies that do 
propagate and contribute to the radiated noise. This mode trapping and up-scattering was first noted in 
the context of tone noise by Topol, Mathews, and Holheubner (Ref. 9), studied and modeled 
extensively by Hanson (Ref. 1), and incorporated in a quasi 3D prediction scheme (TFaNS) by Topol 
(Ref. 5). For certain vane/blade ratios, mode trapping/frequency scattering dominates tone generation 
in fans. A major objective of the present model is to find how important this mode trapping 
phenomenon is for broadband fan noise. 

The next 3 sections of this report are devoted to building a broadband noise prediction scheme using 
the modal scattering kinematics described above. Section 5 develops a notation for modes (via Fourier 
transforms) compatible with the scattering rules of this section. Then Section 6 derives equations for 
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Figure 1 1 . Mode order/frequency plot showing cutoff limits for both the axial flow region (upstream of 
rotor and downstream of stator) and for the swirl region (between rotor and stator). In the swirl region, 
M y = 0.4. In the shaded region, modes can reflect back and forth between rotor and stator but cannot 
propagate beyond the blade rows. 


coupling rotor and stator using a matrix of scattering coefficients and vectors of modal amplitudes for 
the scattered waves. Section 7 relates a standard turbulence spectrum to a source vector, which drives 
the coupled equations, and presents working equations used in the prediction code CupBB. Sample 
calculations are shown in Section 8. 

The prediction scheme is based on the independent mode set property described above so as to 
minimize the required computer resources. Rather than treating broadband prediction as one very 
large coupled problem, it is broken down into many small coupling problems via loops on frequency 
0) o , radial mode order ju, and the mode order shift index q. Each of these smaller analyses finds 
amplitudes of all of the modes and frequencies in a mode subset (as shown, for example, in Figure 9) 
by solving a series of simultaneous equations (derived in Section 6) for blade row scattering. 
Excitation for the system is given by modes of the turbulence input at all the same orders and 
frequencies simultaneously. 
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SECTION 5 

FOURIER TRANSFORMS AND STANDARD WAVESET NOTATION 


The preceding section established the scattering rules for broadband rotor/stator interaction. In this 
section, we develop a notation system for Fourier transforms tailored to those rules and designed to 
facilitate development of the code for noise prediction. 

Recall that the totality of modes that participate in the interaction can be broken down into independent 
sub-sets that scatter only among themselves. And recall from Section 4 that for excitation with a given 
frequency offset co 0 , mode offset q, and radial order /u, the set of permitted frequencies and 
circumferential mode orders is given by 

frequency: (O = 0) o + nB { Q., n = ... — 1, 0, 1, 2,... 

. 1 ) 

mode order : m - q + nB\ - kBj k = ... — 1, 0, 1, 2, ... 

i.e., frequencies are separated by blade passing frequency and circumferential mode orders are 
restricted to the values of m resulting from integer values of the rotor and stator scattering indices, n 
and k as shown in Figure 9. Thus, in the computer code, we can set up loops on aib, p, and q and 
solve a relatively small coupled problem for scattering on n and k for each of these combinations. In 
each case, the input is a vector of turbulence modes for all n and k and the response is found 
simultaneously for the acoustic modes with the same n ' s and k' s. 

The pressure and vortical wave systems are treated separately below. 


Pressure Waves 

Background flow for the unsteady interactions will be divided into 3 distinct regions as shown in 
Figure 3 and counted by the index r : 

r = 1 : upstream of rotor 
r = 2: between rotor and stator 
r = 3: downstream of stator 

Furthermore, in each region, there are 4 distinct wave types (or families) counted by index T : 

7 = 1 : Upstream-going pressure waves 

7=2: Downstream-going pressure waves 

7=3 and 4: Two independent sets of vortical (convected) waves 

These indices are used to express pressure waves (7 = 1 or 2) in region r as a typical Fourier integral 

p'j (x, 0 = £ P T r (X, to) e~ im dco (5.2) 

To coordinate with the frequency scattering rule, we break the frequency axis into segments each 
BPF long 
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(5.3) 


(x,/) 


£ Pf{x,co 0 +nB&)e^<»o + »W'd<o 0 


and then further expand the Fourier transform as follows to coordinate with the radial and 
circumferential scattering rules 


Pl(x,m 0 + nB l Cl) = p 0 ^ £ I P-f^i^+nB^) (5.4) 

fj -0 q-0 k=—°° 

where p 0 is a reference pressure taken to be the sea level, standard day value. Xj is the .v reference 
location for the type T waves in region r. We have an independent scattering problem for each 
combination of co 0 , p ,q. The solution to each of those scattering problems gives Hnqk^o + 
for a range of n's and k' s as shown in the mode/frequency plot of Figure 9. 

We will solve for only a small number of a> 0 's and they will be discretized as coj. For example, if we 
are satisfied to discretize the frequency axis with points at the BPF harmonics and points midway 
between, we would use 

©, =0.5B,£2 

(5.5) 

co 2 = B X Q. 

Along the frequency axis represented below, the simultaneous solution for coi would find the spectral 
levels for all frequencies marked by x and the solution for Oh would find the levels for the 
frequencies marked by •. 

•- — X- — •- — X- — •- — X- — •- — X- — • x •- — X 

0 1 2 3 4 5 frequency in BPF orders 

The velocity components for the 2 pressure wave types (and the vortical waves) can be written in form 
parallel to Eq. 5.3 as follows. 


uj(x,t) = £ t/f (x,^ +«g,ft) e- i( ^ +nB ' C1) ’ dco 0 

n =- <*> 

vf(x,()= £ J 0 S,a Vf(x,< 0 „ +nB l Q) da>„ (5.6) 

^ = -00 

v4 (x, 0 = £ wf (x, 0 O + nB x Cl) e~ i(c °o dc0o 

n=—<x> 

where 
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oo B ~1 oo 

U T J (x, a> 0 + nS,Q) = £ £ £ l/f Mk (<° 0 + «B,Q) ( ' wf )+ '” 01 cos(^) 

ju=0 <7=0 k-~oo 

CO B2 — 1 OO r r 

V{(\, C0 o + nB x £l) = a Q £ £ £ Vf Mk (co 0 + «B,Q) ( '“ v )+,M ^ ] cos(^) 

fj — 0 ^=0 k=—°° 

Wj(\,(0 o + nB\Q) = a 0 £ £ £ Wj Mk (co 0 + «B,D) (A_Ar )+ ”^ ] sin(^) 

a=0 <7=0 £=-00 


(5.7) 


Later, we will need to express the velocity transforms in terms of the velocity perturbations to 
represent turbulence (i.e., to invert Eq. 5.6 for T = 3, 4). These can be found in the usual manner from 
Eqs. 5.6 and 5.7 by exploiting the orthogonality of sinusoidal and exponential functions. For example, 
the transverse velocity transform is 




4n RaJiUSjj 


C°° rxR rh n 

L Lr lo s®* 


-i(k 




cos ( — — ) dx dy dz (5.8) 


where = 1 for ju = 0 and = V> for /u> 0. x =(x, y, z) is the fluid-fixed coordinate. vf is the 
velocity represented in the fluid-fixed frame, i.e., vy(x) = vf (x + Ur,<). Th e transform of the axial 
velocity is the same but with vf replaced by u' T . 


[Note that in this type of waveform analysis, the integral in Eq. 5.8 will not converge if the extent of 
the velocity field is infinite in the .v direction. One procedure for dealing with this is to set the 
perturbation velocity to zero outside of a region given by - Ux < x < Ux where 2x is the time for the 
block of flow to pass the cascade. This guarantees convergence of the integral. When the desired 
energy spectrum is found, it can be divided by 2x to find the power spectrum. Then we can take the 
limit as x approaches infinity.] 

In the following, we define some non-dimensional quantities based on a representative radius R. This 
is a constant for the problem and can be thought of as a scaling parameter that brings the notation for 
our rectilinear representation of the fan closer to notation used in full 3D annular representations. 
Thus, from Eqs. 3.13 and 3.14 


k r = k r R = 

K XTnk ~ K XTnk n 




M x (-<£> + mM ’ y ) + J(-db + mM r y ) 2 - /3~(m 2 + v 2 ) 


(5.9) 


with the upperjower sign going with T = 1,2 and from Eq. 3.7 


AL ~yL = -d ) + M r J r mi +M r m (5.10) 

a r 

with d> = coR/a r =^Y^ + n^B\M’ t where Af/ = ClRja,. is the tip rotational Mach number and M' x and 

My are the axial and swirl Mach numbers. All Mach numbers are based on local sound speed in the 
regions counted by r. 
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In defining the standard wave sets, we deal with defining components and associated components. 
The pressure waves are defined in terms of the pressure amplitude coefficients in Eqs. 5.3 and 5.4. The 
associated acoustic velocity components can be related to the pressure coefficients by inspection of 
Eqs. 3.15: 


Pjfxqki^o + wBj£2) = defining coefficient for pressure waves 
P T M k ( (0 o + nBfl) 


U r TMk ((0 o + nB ',Q)=— ^ 


( -k r ' 

_XTnk 


jlqk (°^o + nBfo) 

^T/uqk^o + nB { Q.) = 


p,.a 0 a r 


_ Po 


A 


Tnk 


Pr a o a r 

Po 


—m 

^■Tnk 


p T M k ( co Q + nBfl) 


(5.11) 


7 = 1,2 


p r a 0 a r 


-mnjH D 


A 


Tnk 


P Tj.lqk( (0 o +nB { Q) 


Thus, if P-fuqkiWo + nB\Cl) were known, the pressure field could be constructed from Eqs. 5.3 and 5.4. 

In fact, in Section 6, the P T M k ( co 0 + nB\Cl) transforms with ranges of the indices r, T, n, k will be 

used as elements of the system state vector similar to the Fourier coefficients in the state vector for the 
tone problem in Eq. 2.5. 


Vortical Waves (7= 3, 4) 

Since vorticity is a 3D vector, it might appear that 3 components could be specified independently. 
However, because the convected field satisfies the continuity equation V • u = 0 , only 2 are 
independent. In fact, for the p = 0 case, which corresponds to 2D flow, there is only one component. 


We deal with the p > 0 cases first and define 2 families as follows 

7=3 family: has v and w components but no u 
7 = 4 family: has u and w components but no v 


The 7=3 family is defined via the transverse velocity component in a form parallel to the pressure 
wave families. Thus, the V^ qk ((D 0 + nB\Q )' s are the defining coefficients (that go into the state 

vector). Again, the W^k^o + ’s are related through the continuity equation, dw/dz = —dv/dy . 

Hence, 


for p > 0, 


p 3 M k(G>o+nB l Cl) = 0 
U^ qk ((0 0 +nB^) = Q 


(co 0 + nBiQ) = Defining component 




^ —im ^ 


pnjH D j 


P 3 nqk (0) o + flB { Q) 


(5.12) 
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The T = 4 family is defined via the axial velocity. Thus, the U^ qk {(0 o + nB x Q.) 's are the defining 

coefficients (that go into the state vector). The W^ qk (0) o + nB^Q.) ’s are related through the continuity 
equation, dwjd: = -du/dx. 


for n> 0, 


Kjdqk ( W o + n B\Q) - 0 

U’*Hqk ( 0J o + nB\Q) = Defining component 


^4 M k(0> o +"B l Ci) = 0 
^4/Uqk («o + nB\Q.) — f 




~ ik ’x4nk 
flK I H D 


\ 

U 4 fiqk nZ?j£2) 

/ 


(5.13) 


Finally, we define the vortical waves for /u = 0. They only exist for the T = 3 type. We still choose to 
define these by the transverse component but now the associated component is the axial velocity via 
du/dx = — Bv/dy and the radial velocity is 0. 


for n = 0, 


p 3nqk(° ) o +nB lQ) = 0 

U’lnqkiMo + nB x Q) = 


r \ 

Vl Mk «o 0 +nB x &) 

( .v3n£ J 


Vlfiqk (° } o + »B\Q.) = Defining component 
w 3 Hqk (Mo + nBlQ) = 0 


(5.14) 


The reader can verify that the wave families defined above satisfy the continuity and momentum 
equations in each of the 3 flow regions and that the boundary conditions of periodicity in (p and flow 
tangency at z = 0 and h are satisfied. 

In the above, wavenumbers for the convected waves are normalized (from Eq. 3.17) as follows 


k x3 nk = Kx4 nk ~ 


G) — mM j 
M r 


(5.15) 


In the next section, we will show how the defining components established above are used in the state 
vector and coupling equations. 
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SECTION 6 

SCATTERING COEFFICIENTS AND THE COUPLING EQUATIONS 

In this section, we develop the rotor/stator system coupling equations; these match the modal output of 
the stator to the modal input of the rotor and vice versa and show how the turbulence enters as a 
source, i.e., excitation, for the system. This is all in terms of Fourier transforms of the acoustic and 
vortical perturbations - as if they were known. Since, in fact, they are knowable only in a statistical 
sense, both the turbulence and resulting noise must be treated in terms of expected values of 
magnitudes of the Fourier transforms. That subject is treated in Section 7 where we show how to 
compute the power spectral density of the sound from the turbulence spectrum. 

In the preceding section, we developed expressions for all permitted perturbations in the 3 regions 
under consideration: 

Region r - 1 - upstream of rotor 
Region r = 2 - between rotor and stator 
Region r = 3 - downstream of stator 

The waves are of 4 types: 

Type T = 1 - upstream going pressure waves 
Type T = 2 - downstream going pressure waves 

Type T = 3 - vortical waves (first independent set based on v — velocity component) 

Type T = 4 - vortical waves (second independent set based on a - velocity component) 

Perturbations were expressed in Fourier expansions. If the Fourier coefficients were known, the flow 
would be completely defined from Eqs. 5.3 to 5.14. The coefficients defining the waves are of the 

form Pf^ q k (&> 0 + nB\Q .) . We form a state vector of these coefficients and compress the notation for 

convenience. Thus, for each region r , 


A[(n, k) = P{ Mk (Q ) 0 + «5,£2) 

Aq_ («, k) = P{ M k (0J o + nB l^ 

A$(n,k)=V 3 r Stlqk {(D 0 +nBp.) 

A 4 (n, k) = U\ Mk (co 0 + /jfi,Q) 

The ordered array of the (n, k) ’s is the state vector 

A. The n’s run from 1 to a user-specified upper 
frequency limit and k" s run over ranges corresponding 
roughly to cuton for the acoustic waves and to the 
turbulence spectrum for the vortical waves. The /J's, 
q' s, and (0 0 's are not needed in the state vector 
notation because they are constant parameters for 
each coupling sub-problem. 



( 6 . 1 ) 




To help understand the structure of the coupling 
equations, consider the scattering of modes by 


ROTOR 

Figure 12. Sketch to explain rotor scattering 
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the rotor shown in Figure 12. We focus on the output vortical waves (Type 3) in region 2 (to the right 
of the rotor) and what contributes to them. These can be scattered into from each of the 4 input wave 
types and that action is denoted by scattering coefficients as follows. 


A^{n',k) = X [S$}(n',k;n,k)/f(n,k) + sl\{n\k-n,k)A\{n,k) 

n 

+S33 ( n ,k\n , k)A\(n,k) +534 (n',k’,n,k)A^(n,k) + S3 (n', 


( 6 . 2 ) 


This meaning of this equation is that the Type 3 vortical waves in region 2 come from scattering of the 
up-going pressure waves (T = 1) in region 2 and scattering of the Types 2, 3, and 4 waves in region 1, 
The summation is on the n index since we found in Section 4 that rotors scatter on n and not on k. 

We have also added S3 (n',k) . These are source waves that are to be prescribed as excitation for the 
system. Notation for the scattering coefficients is 



Figure 13. System coupling matrix 
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f . 

Sj^r{n',k'\n,k) with the unprimed indices denoting input modes scattered from and the primed 
indices denoting output modes scattered into. 

Equation 6-2 can be generalized to treat the entire inlet/rotor/stator/nozzle problem as shown in 
Figure 13. Equation 6.2 can be seen as the seventh line of the matrix equation. All of the other 
scattering coefficients are shown as well. Dots indicate empty elements, or no scattering. The next 
sections show how the stator and rotor scattering coefficients are computed. (We ignore scattering 
from inlet and nozzle in this report but expect to address it in the future.) 


Computation of Stator Scattering Coefficients 

The stator acoustic element consists of the stator vanes themselves and an actuator disk which 
represents turning of the mean flow. The jump in mean flow properties across the disk determines 
jumps in perturbation quantities via conservation of mass and momentum. In this section we first treat 
scattering by the vanes and then later show how the actuator disk effect is added. 

To compute scattering by the stator vanes, note that the connection between input modes and output 
modes is through the upwash at the stator. Thus, we need to Find the upwash associated with each of 
the 4 input wave types. Then, we need to find the amplitudes of each of the 4 response modes to unit 
upwash. The required scattering coefficients are formed from the 16 combinations of input and output. 
The upwashes associated with various wave types follow from the wave kinematics of Section 5 and 
are shown on the right in Table 1 under “Excitation.” The output mode amplitudes due to unit upwash 
are found using Glegg’s Wiener-Hopf cascade method (as explained later in this section) and are 

r __ 

shown on the left in the table under “Response.” Notation Sj^(n,k';n,k) is used for scattering by the 
vanes. Its meaning is given by 

Aq->(n,k ) Sj>f(n,k \n,k)Aq-{n,k) (6.3) 

When the actuator disk is added, the coefficients will be modified and placed in the scattering matrix 
with the same notation. For example, scattering of vortical waves into downstream pressure waves is 
given by 

A 2 ( >i,k )< — S 2 ^(n,k ;n,k)A i (n,k) (6.4) 


The prime on the superscript on the LHS indicates that the coefficient does not yet apply to Region 3; 
this will be dropped when the actuator disk is included. This coefficient for scattering of the type 3 
vortical waves (with fj. = 0) into down-going pressure waves can be constructed from elements of 
Table 1, 


S 2 }(n,k';n,k) = 


1 Ain c 2 P 2 A IndknkLPziM' \ m 


82 


R 


nk' 


Y W 

^ ^.v3 nk 


-sin0 2 +cos0 2 


(6.5) 


Derivations of the excitation and response factors in Table 1 are discussed in the next 2 sections. 
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TABLE 1 

FACTORS FOR STATOR SCATTERING COEFFICIENTS 


T 

Response 

T 

Excitation 

1 

a 2 , _ — 1- 4 / /r <T 2 Pi ^Ink'Clnk' D \nkk' w nk 

1 

nVi* ^2 ^.vi«A sin 0 2 ~ w cos 0 2 ^ 3 ' , , x 

A l\ n ’ k )~ - . „ 

82 R>ik ' a o 

«. " ^ Af, t M 

2 

a 3 ',„ _ +l'4/^r 2 P 2 A 2 n k’ Clnk’ D 2 nkk' w nk 

2 

w nk _ (? 2 KA sin 0 2 -WCOS 0 2 2 M 

a 2 (n,A ) - _ 

'V 

^ L4/ 2 A 2 „ a . 

3 

all n 

,3'/ > —2ki 

^3 («.*)=_ 

g 2 cos &2 

^ 'f £ - sin 0 2 + Cc cos ^2 1 

k ? f 2 +«„V+v 2 J 

K nk 

a o 

3 

t^O 

( \ 

— sin(? 2 +cos0 2 A 2 {n,k) 

a ° k ^.\3nk j 

n>o 

^lL = (cosd 2 )A 2 (n,k) 

a 0 

4 

u>o 

„3\ , ,, ~ 2 Ki 

A 4 (n 9 k ) - _ 

cos 02 

cos0 2 ~Cc sin ^2 " 
fc+4'+v 2 , 

y ^nk 

“ nk 

a o 

4 

/J > 0 

= (- sm 9 2 )A%(n,k) 

a o 


Above, is the potential jump across the wake from Glegg’s theory (see Appendix A). Subscripts nkk' 

indicate input waves with interblade phase 


a = • 


2jzm 

~b7 


and output waves with interblade phase 


m = q + nBi — kB 2 


Also, 


/ 

a 


2k m' 


m' = q + nBj — k B 2 


A = — d> + M x k 2 T ' + M y m' 


and in 


Ct nk ' = ~V ± (ft 2 m' + 0)M y )C0S 0 2 ] 

Px 


the upper, lower sign goes with T' = 1,2. K„k is Glegg’s vorticity factor (K(y c ) in Appendix B), which depends 
only on n and k (noton k'). The wavenumbers required for vortical response are 


f c =d)/M 


«„*' = sinS 2 + m'cosSj 



C 2 =-tf?+v 2 ) 
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Excitation Factors 

Here, as an example, we show how the excitation factor for down-going pressure waves at the stator is 
derived; this is the second item in the right column of Table 1. The remaining excitation factors can be 
derived in similar fashion. 


The upwash at the stator due to the down-going pressure waves in region 2 is given by 


m’(x, t) = — 1*2 (x, t) sin 0 2 + v ’2 (x, t) cos d 2 


( 6 . 6 ) 


By using Eqs. 5.6 and 5.7, it can be shown that the coefficient of the upwash is 

w nk = a o 2nqk (<9o + nB^Q) sin 0 2 + V 2/Uqk (0) o + nB x Q) cos 0 2 J (6.7) 

The U and V coefficients are related by Eqs. 5.9 to the P coefficient that defines the pressure wave. 
When these are substituted into Eq. 6.7, the result is 


w nk _ Po 


a o Pl a 2 a o 


, A 2 


khnk sin 0 2 - m cos 0 2 


A 2 nk 


l/uqk + nBlQ) 


( 6 . 8 ) 


The factor before the square brackets can be modified using p 2 a 2 = yp 2 where y is the ratio of 
specific heats (y= 1.4 for air). Using this and replacing P 2klqk (C0 o + nB i Q.) with A 2 (n,k ) lead to 


m '»A- _ ^2 

1-4 P 2 


sin6> 2 — wcos0 2 


A 2 (n,k) 


(6.9) 


which is the second item on the right in Table 1, as desired. A ? and P 2 are speed of sound and pressure 
in region 2 normalized to sea level, standard day conditions. 

Response Factors 

These are the factors in the left column of Table 1 and are computed using Glegg’s Wiener-Hopf 
cascade method. As an example, we derive the modal coefficient for upgoing waves from the stator 
due to unit upwash (the first item in column 1 of Table 1). 

Appendix A adapts Glegg’s equations to duct coordinates and to our modal indices n, k , and p. 
Equation A-43 gives the pressure wave scattered into mode n,k' due to upwash in mode rt,k as 


Pnkk’ 


-in p 2 a 2 c 2 
8 2 


■ w nk 


A 2 

A 1 nk' 


Clk' 


D 


Inkk' 


R 


i[k 2 xlnk -x+m'<t>-(Ot] Mnz 

e cos(— — ) 


nk' 


( 6 . 10 ) 


where we have replaced Glegg’s ± superscript notation indicating up/down going waves with a 
subscript 1, since we are addressing upstream going waves. Upon comparison with Eq. 5.4, we can 
write for the modal coefficient of the upgoing wave, at frequency a>-0) o + nB^Q. 
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( 6 . 11 ) 


P\ M k((0 o + nB^) = 


_ rn P 2 a 2 C 2 h\ nk ’ C\ nk ' D[ n kk’ 


Po8 2 


R 


nk' 


w nk 


Note that the phase in Eq. 6.1 1 is based on an .v-origin at the leading edge. In code CupBB, origins are 
shifted to the appropriate reference planes after coupling to the actuator disks. 

In Eq. 6.1 1, we again apply p 2 a\ — 1.4/>2 an d replace P^picOo + n B\Q) with Aj 2 ( n,k ') , leading to 

A~(n k') — EZ ^2 ^lnk' C\nk' P^lnkk' w nk ^ 

g 2 A 2 R nk ' a 0 


This is the response factor giving the modal coefficients for up-going waves from the stator due to unit 
upwash. It is the first item in the left column of Table 1, which was to be found. The factor relating to 
downstream pressure wave response follows immediately. The response of the vortical waves can be 
found from Equations in Appendix B. 


Rotor Scattering Coefficients 

Derivation of the rotor scattering coefficients proceeds in the same manner as that for the stator. The 
16 groups of coefficients are formed from factors in Table 2. The 4 inputs in the right column drive the 

4 responses in the left column. The resulting coefficients are denoted by Sj^{n',k\n,k) . They are then 
modified by the effect of the rotor actuator disk and placed in the system scattering matrix S. 
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TABLE 2 

FACTORS FOR ROTOR SCATTERING COEFFICIENTS 


T’ 


Response 


Excitation 


_ — 1-4/ 71" q P I A in’k Cin'k D\„kn’ w nk 
Aj {n ,A) — — 


Si At 


R, 


w nk _ A 2 sin - m cos ^ (/|> , ;) 




1.4P 2 


MnA 


A 2 (n / A') pi A 2 „ ' k Cln'k R*ln kn' ^nk 


Si 


R 


w nk _ A 2 ^ sin6>i - WCOS 01 


n'k 


a 0 1-4P- 


aLa- 


all t< 


U=0 


a|(«U)=— 




£l cos0j 


^ 7, A sin cos N 

?c + <& +V 2 


H’ i 


nk __ 


AT; 


w. 




ajA 




m 


\y3/jA 


sin 6^ +COS0J 


A 3 (/i,&) 


ff/, 


i />0 


_ ( ^ \ ><1/ 


• = (cosfyMj («,£) 


jU>0 

,1 2 / ' / \ 2tt/ 
44 (« ,£)=— 

& 1 cos 0j 


iU>0 


^ A 


<V* ft- cos 0i -C 2 sin0| 

£+<&+1? 2 




'nk 


w nk _ , a \ a l' / 


• = (-sin#!)/^ (n,k) 


Above, Drnkn' is the potential jump across the wake from Glegg’s theory (see Appendix A). Subscripts nkn ' 
indicate input waves with interblade phase 


<7 = ■ 


2nm 

~b7 


and output waves with interblade phase 


m - q + nB i - kBj 


, 2nni 
O = 

Bo 


m — q + n'B\ — AB 2 


and 


A r'n'k =-<*>+ M x P rT ' n ' k + M yin' 

£ TVk sin0 i ±(Pfm'+<a r (M y -M T ))cos0, J 

Px 


where /3 2 = l-M? and M, is relative Mach number at the rotor and the upper,lower sign goes with the 
T=l,2wave. K nk is Glegg’s vorticity factor (K(y c ) in Appendix B), which depends only on n and k (noton 
n'). And the wavenumbers required for vortical response are 


f c = <£>,. /M l a n ' k = -k 2 3n > k sin d { + m cos 9 { v= ^ K / Hd £c = ~(?c + v ~ 2 ) 
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Combining Blade and Vane Elements with Their Actuator Disks 

The coupling scheme shown in Figure 13 treats 4 acoustic elements: Inlet, Rotor, Stator, and Nozzle. 
Only the rotor and stator are addressed in this report; inlet and nozzle reflection and transmission must 
be treated in a separate study. Figure 14 shows the mean flow entering the fan stage axially, turning at 
the rotor, and straightening back to axial at the stator. To deal with this turning, actuator disks are 
combined with the rotor and stator blade rows. Appendix C derives the actuator disk equations based 
on linearized equations for conservation of mass and momentum; jumps in background flow properties 
cause jumps in the wave modal amplitudes. 

Code CupBB models the background flow via ID isentropic (lossless) flow equations and the standard 
equation for action of the rotor 

C„ 

A(tangential velocity) = x A(total temperature) (6. 1 3) 

wheel speed 

Ambient conditions correspond to a standard day at sea level. The user inputs inlet axial Mach 
number, rotor rotational Mach number, and stage pressure ratio; the code calculates the remaining 
mean flow quantities based on isentropic flow and Eq. 6.13 assuming axial flow out of the stator. 

This section shows how the actuator disk scattering factors are combined with the scattering factors 
just found for the blade rows to produce a combined rotor/actuator disk element and a combined 
stator/actuator disk element. These are the elements that are used in the coupling matrix of Figure 13. 

Figure 14 shows the rotor with its actuator disk, the stator with its actuator disk, and all of the wave 
families included in the coupling scheme. As described earlier, the A’s are elements of the state vector 
(to be found) and the B' s are elements of the source vector (to be prescribed from a turbulence model). 
The A’s and B's are further subscripted on the scattering indices n and k. The derivation below is 
in the context of the stator but the method of combining a blade row with its actuator disk is the same 
for the rotor as for the stator. 

The lower part of Figure 14 shows the stator element broken down into the vane row and the actuator 
disk. Scattering equations for the vane row are 

A? = S°2 4 + Sf" Af + S™A% + Sff Af + 

4 = S&ti + S%A$ + S%4 + S%4 + B b 2 (6 14) 

4 = 4 i a i + sf? 4 + sff 4 + s* 1 4 + 4 

4 = s&4 + sjgAg + s b 44 + s!fi4 + 4 

and for the actuator disk are 
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(6.15) 


A? = Sg4 + S$A$ + Sft4 + sff Af 
Aj = ^22^2 + ^23 A3 + S74 A4 + Sfl Af 
A 3 = $32 A 2 + $33 4 + S34A4 + ^3 ! A \ 

4 = s $4 + 543 A3 + S$a£ + S44 Af 

These can be placed into a matrix system similar to that in Figure 13 as follows 



(6.16) 


With the block matrix notation indicated by the shaded boxes Eq. 6.16 reduces to 


A 0 = $AA A ° + $ AB^ 

+B a 


^ ~ $ BA A? + $BB^ b + $BC^ C 

+B b 

(6.17) 

J? 

+ 

8 

to 

11 

+B C 



The middle line of Equation 6.17 can be written 

(1 - S BB )A b = S BA A a + S BC A c + B b (6.18) 
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and solved for A b via 


A b = E S BA A a + E S BC A c +EB b 

where 

e = (i-s bb t 1 


(6.19) 

( 6 . 20 ) 


Substitution of Equation 6.19 into the first and last lines of Equation 6.17 eliminates the scattered 
waves in the middle ( b ) region 


or 


A* = + S AB E S BA )A a + E S BC A c + B a + S AB E B b 

A c = S CB E ^BA^ a + (^CC + Scb E S B c)A c + S CB E B b + B c 


( 6 . 21 ) 


V 


s aa +s ab E S BA j S AB E S BC 

X 

A a 


B a + S AB 

EB h 

A c 


Scb E S ba j S cc + S CB E S BC 

A f 

+ 

S CB EB b 

+ B C 


( 6 . 22 ) 


In the current application, we will not prescribe any source waves in the “fo” region so that finally 


A a 


5 A4 ES ba 

Sab E S bc 


~A a ~ 


B a 

— — 

— 

— — 

X 

— — 

+ 

— 

A c 


Scb E S ba 

s cc + Scb e Sbc 


A c 


B c 


(6.23) 


This is the desired form for combining the stator vanes and the stator actuator disk into a single 
acoustic element. The scattered waves in the region between the vanes and the disk have been 
eliminated from the equations so that the input/output properties of the stator element are completely 
specified by the A a and A c state vector components. To apply this scheme in the computer code 
CupBB, we first compute the scattering coefficients for the vanes with uniform background flow 
corresponding to region 2 (between blade rows) and place these into the S array as in the example of 
Eq. 6.5. Then we compute the scattering coefficients for the actuator disk according to theory in 
Appendix C and place them in array D. Next, elements of S are placed into block matrices S 'aa, S AB , 
S B a . S B b and elements of D are placed into S BB < Sbc< $CB’ See per Eq* 6-16. Equation 6-23 is applied 
to find an array to replace the old array S. Finally, sections of S are placed into the stator location of 
the system scattering matrix of Figure 13. The same procedure is applied to the rotor blades and their 
actuator disk to complete the system scattering matrix. 
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SECTION 7 

EQUATIONS FOR SOUND POWER SPECTRA 


This section develops equations required to compute spectra of sound power propagating upstream and 
downstream from the coupled blade rows. The development is in 5 stages: 

1 . Derive equations for sound power in terms of the Fourier transforms defining the standard 
wave sets in Section 5. 

2. Derive equations to compute those Fourier transforms in terms of the inverted coupling matrix 
and the source vector elements. 

3. Develop equations for the source vector expected values in terms of 3D turbulence spectra. 

4. Adapt a standard turbulence spectrum (the Liepmann spectrum) to the ducted flow 
environment. 

5. Assemble the working equations for programming. 

Sound Power Equations in Terms of Standard Wave Set Coefficients 

Power associated with acoustic perturbations is evaluated according to the energy flux vector 
given by Goldstein (Ref. 10). 


1 = 


■ + u-U 


Pr 


(p,.u + pU) 


(7.1) 


This is the instantaneous energy per unit area per unit time in the direction given by I. The perturbation 
pressure, density, and velocity are p, p, and u. p ; and U are the density and velocity of the background 
flow in Region /*. The energy flux vector is space and time dependent. We are interested in its 

component in the ,v direction. Forming I x = I i x and applying p = p/a r lead to 


, M 

I x =(1 + M 2 x )pu + £■ 

p,.a r 


2 2 

p + p r a,.M x u + p r a r M yUv + M X M y pv 


(7.2) 


for the power flux in the x direction. M x and M y are Mach number components in the axial and 
tangential directions relating to the background velocity U = ( U , V, 0). 

Recall from Section 5 that the expressions for perturbation acoustic quantities are 
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pf(x,/) 

u’ T (x,t) 

°° - t? O 00 B 2~\ 00 

-ft IT I 2 I 

n ~ — oo u =0 q = 0 k=—°° 

i 

k.XTnk 

Pr a r^Tnk 

vf(x,r) 


— m 

L J 


_ Pr a r^Tnk _ 


P Tjlqk (®o + n B\Q) 


Xe M.\Tnk (-v-.v 7 r )+m<t>-(a> 0 +nB l Q.)t] 


COS 


r jJL TTZ 
x h 


\ 

) 


d(0 o 


(7.3) 


Before substituting these into Eq. 7.2, note that, in our convention, the pressure and velocity are pure 
real by virtue of using 2-sided series (i.e., imaginary parts in the upper and lower halves of the n series 
cancel). Thus ,p 2 =pp*, and we use pp , for convenience in the manipulations that follow. 


Lx = 


_ Po 


oa 1 1 JJ III III 

r n n p q k p q k 


(1 + M~) 




+ M x + M x 


( -k r V 

'hank 


l Tnk 


-k ,r ' ■ 

K XTn k 

A rr 

Tnk' 


+ M y 


-k 


XTnk 


V -m' ' 


A", 


a;: v 


+ M X M y 


( -m X 

^Tnk jj 


^TjJqk^O + n B\Q) Pt jj'qYit&o + H ^ 1 ^) 


xe 


y Tnk Jy Tnk ) 

il( k ATnk )+(m-m')<l>-[(cQ 0 +nB x £l)-((d' 0 +n%£l)]t ] Q f P Kz 


l * 


:cos 


T 


- \ 

) 


d(o 0 dco' 0 


(7.4) 


We need to average this over time and then over the space variables r and <p. For the time average, 
recall that the Fourier transforms apply to a block of flow passing through the cascades in the period 
-r < t < T. Outside of that time range, the perturbations are temporarily set to 0. We still integrate 
over all time since it produces a Dirac delta function 



(7.5) 


This gives the energy due to that block of flow. Then we divide by 2 r to get the average energy/unit 
time, i.e., power. Similarly, the average over r requires 


if 

h J o 




(7.6) 


where £ u - 1 for p = 0 and = Vi for p > 0. Also, the average over 0 yields 5 qq ' <5«-. These operations 
eliminate all the cross terms in Eq. 7.4 with the result 


1 1 III 

Pr a r n \x q k 


K£, 


P 1 

r Tmk 


(0) o + «B!0)| 2 


dco„ 


(7-7) 
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where 


F nk ={\ + M x ) 


r -k r ' 

K .\Tnk 


V IV Tnk ) 


+ M x + M x 


k'xTnk 

^■Tnk 


+ My 


K XTnk 


m 




'■Tnk 


A 


Tnk 


+ M X M y 


( \ 

—m 


A l Tnk 


By applying definitions of the wavenumbers in Eqs. 5.9 and 5.10, F nk can be reduced to 

+P x (b>[E 


F„k = 


(-&> + mM y ) + M x Je] 


where E is the cuton discriminator 


(7.8) 


(7.9) 


E = (-cb + mM v ) z - p~ (m 2 + v z ) 


(7.10) 


This is positive for this application, since only cut on waves carry power. The + sign in the numerator 
of Eq. 7.10 indicates that energy upstream travels in the negative x direction and downstream travels 
in the positive .v direction. This sign is dropped henceforth. 

Now, since the expression for I x is an integral over frequency, we identify the integrand as the power 
spectral density of acoustic intensity. (We still have to multiply by duct area and bandwidth to get 
sound power.) We define the area averaged intensity spectra upstream and downstream as 


U l n q (o ) 0 + nB^) = 
Yl] lq {(0 0 + nB x a) = 



P\ a \ 

Po 


P 3«3 



7ZE, 


1 Mt 


{(D 0 + 





( (O a + nB { a) 


(7.11) 


for each independent mode sub-set defined by 0) o , q, and fd. Next, we will find the required expressions 


for 


i 2 , 2 

Py k (0) o + nB { Q) and P 2/Jqk ( co 0 + nB { Q) in terms of the source vector elements. 


Fourier Coefficients in Terms of the Coupling Matrix and Source 

The terms P* k ((0 o + nB l Q) and P 2(iqk (co 0 + nBjQ) are elements of the state vector developed in 
Section 6. Specifically, 


P y \ qk (0) o + nB { Q.) = a\ ( n,k ) 

and (7.12) 

p 2 i 9 A-(®o +nB^) = Az{n,k) 
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are the coefficients for upstream-going pressure waves in region 1 and downstream-going pressure 
wave in region 3. The coupling equation was shown to be 

A = S A + B (7.13) 

where A is the source vector, S is the scattering matrix, and B is the source vector. This must be 
inverted to find A as a function of B 


A = (1-S) -1 B 


(7.14) 


In this section, we write out the equations for the desired elements of the state vector in terms of the 
elements of B that relate to turbulent flow into the rotor and stator. In the following section, we will 
express those elements from B in terms of the turbulence spectrum. 


We treat turbulence entering the rotor and stator as the source. The relevant elements in the source 
vector are 


B\(n,k) - \% qk (co 0 + nB\Q.) 
b\ ( n , k) = U 4 n qk (0) o + nB\Q.) 


which are coefficients of the 2 vortical wave types (3 and 4) in region 1 (entering the rotor) and 

B$ (»> k) = V^ qk (co 0 + «S]Q) ^ 

Z?4 (ft. A') = U Zjiqk (®o + nBfil) 


which are the same for vortical flow entering the stator from region 2. (For radial mode index /u = 0, 
the Type 4 waves are absent.) We denote the matrix of the inverted system by C where 

C = (l-S ) _1 (7.17) 


In terms of the elements of C, the coefficients for the upstream-going pressure waves become 


p \ M k (<°o +nB l Q) = II {C(l, 1, 3, k‘ ')B\(n',k') + C(l,l,n,k-,lA,n',k')B\ {n'X) 


/ i / 

n k 


+ C(lXn,k;2,3,n,k')Bz(n\k') + C(l,l,n,k;2,4,n,k')BZ(n,k')} 


(7.18) 


where the B 4 terms are absent in the case of /t = 0. We will not bother to write out the corresponding 

l l2 

form for downstream going waves until the end of this section. We can not know P JMk (0) o + nB l Q.) 
in general but we can deal with its expected value (or ensemble average) / Pj Mk (0) o + nB i Q) 


Squaring the magnitude of Eq. 7.18 results in many cross terms which we will justify eliminating. The 
cross terms between regions 1 and 2 vanish because we assume the turbulence entering the rotor to be 


NASA/CR— 2001-211136 


44 



independent of that entering the stator. Furthermore, cross terms between wave types 3 and 4 vanish 
since they were defined as independent waves sets. This leaves 


<'A +«b,£2) 

Illl 

C(\Xn,k\\AyX)C\lXn,k\\AyX)(B\{n,k')BX{n'X) 
C(\Xn,k\2,iyX)C*{\Xn,k\2,iyX)(B]{n\k')Bl\n,k”) 

+ CXlnXXA,n'X')C\\XnXXA,n,k r ){Bl(n',k')Bf (h',A')}] 


+ i 
+ i 


(7.19) 


In the following section, we will see that the n'n" and kTc" cross terms also vanish. 

Source Elements in Terms of Turbulence Spectra 

We will treat one of the source expectations on the right of Eq. 7.19 and then write down the rest by 
inspection. Including the 7t8|j/T from Eq. 7.7, we use A to denote 


A - ■ 


KE 


/U 


(7.20) 


T \ ' ' It 

and we proceed by inserting the inverse transform of the F’s given in Eq. 5.8. 


A = 


- 12 


An" Ra 0 hU £ ^ 




(7.21) 


xcos { BEjl cos tULl. dx dy dz dx' dy' dz' 

{ * J l h 


and change variables via 
so that 


x' = x + s 


X = X + s x 

y=y + s y 

z = r + 5- 


(7.22) 

(7.23) 


Then, 
A = - 


n 


JJJJJJ (v3(x)v'5(x + s) 


{An 2 Ra 0 hU) 2 t£^ 

xcosfMiTcosf^lcos 

1*1 1*1 


^x3 n w k*}* ^.v3 n*k wS x+{W m ^R m 


'r 


f fj.ns. ^ 




dx dy dz ds x ds y ds. 

(7.24) 
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The sine term will vanish in the z integration. In flow problems with no boundaries, the standard step 
at this point would be to recognize the expectation of the v product as the correlation function 

(v](x)v](x + s)^ = .K* v (s) and claim it to be a function only of the separation vector s (and not of the 

position vector x). In bounded flows this is an approximation that we assume to be acceptable since 
the length scale turns out to be small (typically 2% to 4% of fan radius). 

We still have to deal with convergence of the x integral. Recall that we are considering only a block 
of flow that is finite in the x direction and then taking the limit as the block size increases to infinity. 
For this, we employ a “top hat” function J(x) which = 1 for -Ut < x < Ur and 0 otherwise. Then 
the correlation is written 

(v 3 (x)v](* + s>) = R l vv (s) x J (a*) (7.25) 

Use of this function permits us to separate the separation and position integrals as follows 


A = - 


7T 


-JJJ cosfek ds y ds. 


(47l-Ra 0 hUf£ v 

A '„R l ,UT . 

JO i-nR f j-t/r { h ) 


(7.26) 


In the first line of Eq. 7.26, we can replace the cosine with an exponential because the correlation 
function is even. In the second line, the z integral yields h D e^ so that 


A = - 


m " HR ,, 

K CCC 1 d .v3 n w k w ^x^ nT^v"* 

7737 /If ^vvO S)e R - h ds x ds y ds : 


(4n RaJJyhx 


(7.27) 


[ UT Uxin'k ■ )X+(m-m')± ] ^ 

J-xR J-Ur ■ y 


The y integral produces a Kroneker delta 


f e R dy = 2nRS n 

J-nR * ” 


(7.28) 


so that the x integral becomes, after applying Eq. 3.17 for k x , 




1 rU X -i(n'-n')- 

- f e u d.x = 2U5„'„’ 

x J-U T 


'n n 


(7.29) 


in the limit as r approaches infinity. With these results, A further reduces to 
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A = 


8 


'n n 


k'k’ 


An-Raluh 


I 


^ m un _ 

i l(kx3n’k' s x^ ~n~^y "* ]^ S - 

Ry V (s)c R h ds v ds v ds. 


J y 


(7.30) 


If the fluid medium were unbounded, this integral could be recognized as the 3D turbulence spectrum 
from the standard form 


J /?*,. (s) c' ka ds = {In ) 3 <(>I y (k) 


(7.31) 


We assume that this is an adequate approximation for the bounded case under consideration, so that 

ne. 




(7.32) 


where the wavenumber components are 

k,=k' x: 

a 2~ /r 

u _ Idn/ 

^3 - //, 

Finally, we normalize the turbulence spectra according to 


^ v (k) = /? 3 W 1 2 o!„(K) 
4(k) = /?V 1 2 oL < (K) 


(7.33) 


(7.34) 


for the 2 turbulence components entering the rotor from region 1, and 

<^(k) = flV 2 2 0 2 ,(K) 

*L< k) = /?V 2 <Dj M (K) 


(7.35) 


for the 2 components entering the stator from region 2. Also, we have normalized the wavenumbers 
by R so that the axial wavenumbers 




= 


k x3nk' R = k x3nk' 


(7.36) 


for regions 1 and 2, respectively, and the circumferential and spanwise wavenumbers 


K 2 = w 7 = q + n'By — k'B 2 



(7.37) 


Now we can write down the 4 required source terms. Since we only need the n' = n" and k' = k" 
terms, we eliminate the Kroneker deltas. Furthermore, we drop the primes completely for reasons to 
be seen below. For region 1, 
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and for region 2, 


*g/i 

T 

r 



(«,*) 



|2 


OU) 


— o!,.(K) 

// D COS0J a Q 


2kM 1 A 1 R 
H d cos0j a Q 


4>i„(K) 


X 



IkMiAi R_^l (K) 
H D cos6 2 a 0 vv ; 


re , u 

X 


\Ulnq tk( n ’k) 


2kM 2 A 2 R 
H d cosd 2 a 0 


®J,( K) 


(7.38) 


(7.39) 


where H D = /;//?, Mi =W ;! a i and M 2 =W 2 /a: are the Mach numbers in regions 1 and 2, A / = a t !a 0 , 
and /Ij = a 2 la 0 . 


Choice of Turbulence Spectrum 

Since we do not have any algebraic form for turbulence spectra in ducted flow, we will simply adapt 
the Liepmann spectrum (Ref. 1 1) that has been used with some success in earlier studies. Unpublished 
analyses by Pratt & Whitney show that this spectrum matches wake data remarkably well. The forms 
for the Liepmann spectrum are 


and 


0««(k) = 


0vv( k ) = 


2 » 2 ( 5 k 2 + A-3 2 

* 2 [ l + f 2 (* f +*|+* 3 2 )] 3 

2v 2 ( 5 7-f + kj 

[l + f 2 (A 2 + k 2 + ) J 3 


(7.40) 


(7.41) 


These are such that, when integrated over the 
3 wavenumber components from -°o to -h», the 

results are u 2 and v 2 . However, in our application, 
we use a sum over discrete wavenumbers (modes) in 
the spanwise direction. Since the spanwise 
wavenumbers take on only positive values, we 
must double the sum to maintain the correct integral 
property. A more accurate statement of this can be 
made by inspecting Figure 15. The full area would 
be approximated by summing all of the cells. 

However, with a single sided series, we use just the 
shaded portion. For this, we take Vi the ji = 0 term 
and the full values for the fj > 0 terms. But to double this, we will actually multiply the fJ=0 term 
by 1 and the jj > 0 terms by 2 and sum to obtain the desired result. 



NAS A/CR— 2001-211136 


48 




Working Forms for Sound Power Spectra 

We are now prepared to assemble the results from the above sections to obtain the forms for 
programming. Sound power spectrum levels are obtained from Eq. 7.1 1 by multiplying by duct cross 

section area 2 nRh, multiplying by bandwidth A co 0 = 2nBW , and dividing by the power reference: 
10" 12 watts. Furthermore, since we use only positive frequencies for calculations, we multiply by 2 to 
account for the lower half of the spectrum. In the final results presented below, the use of primes is 
changed in the notation to bring it into conformity with Section 6. Now, a primed mode order is one 
scattered into and an unprimed order is one scattered from. 

The result for upstream sound power is 


PWR.l a ((0 o + = 5.8838x10 


14 


'fiq 

i2 


..... 

r BW^ 

l in - J 

(ft J 


{ I [l C(l,l,^A';l,3,H,A)| 2 oJ, v (K) + |C(l,l,n , ,A';l,4,n,A)| 2 <J>) (W (K)1 

l P l cos6 l k . n k L 

^rX^VX I [| C(l, 1, A'; 2, 3, n, A )| 2 O 2 , (K) + |C(1, 1 , n , A'; 2,4, n, A )| 2 O 2 , , (K )1 ] 
P l cosd 2 k , „ k L J J 


(7.42) 


Power level in dB is lOlogio of this quantity. A, denotes ambient speed of sound normalized by speed 
of sound for a sea level standard day. P, is similarly defined for pressure. 

The corresponding form for downstream sound power follows immediately: 


PWR 3 m ( co 0 + riBfo) = 5.8838 x 10 


14 


r R f 

r BW ' 

, in - ) 

Hz 

/ 


x( X I” |C (3, 2 , n , A'; 1, 3, n, A )| 2 O J. v , (K) + |C (3, 2, ri , A'; 1, 4, «, A)| 2 (K ) 

L P 3 cosd l k , n k L 

+ X |C (3, 2, ri, k'; 2,3, n,k )| 2 (K) + |C (3, 2, ri, k'; 2, 4, n, A )| 2 (K)l } 

P 3 cos6 2 f ~ T L JJ 


(7.43) 


The working forms for the turbulence spectra are 


4>v,(K) = ^ 

<(K) = r„ 


v 2 2 L 5 

Kf + K 3 2 



W 2 n 2 
V / 

( 1 S A 

u 2 2 L 5 

[l + L 2 (tf?+tff +fff)] 
k\ + K% 

3 

y •* 

W 2 K 2 
\ y 

[l + L 2 (K? + K%+K$) 

3 

- 


^0) _ £1 
K \ - k x3nk 


K 2 = m = q + nB[ - kB 2 



(7.44) 
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for flow entering the rotor from region 1 . Note that the intensity and scale for region 1 are to be used. 
£fi= 1 for jj = 0 and =2 for jj> 0. L is turbulence integral length scale normalized by duct effective 
radius R , i.e., L = UR. 


For turbulence entering the stator from region 2, 


4>n'(K)=^ 

0?„,(K) = ^ 


v 2 2 U 

Kx + Kl | 

W 2 K 1 
\ / 

u 1 2 L 5 

[l + L 2 (Kf +K%+K$) 

L kI+kI . _ 

3 

W 2 K 2 
\ / 

1 + L 2 (K? +K% +K$) 

J 


( 2 ) _ £2 


K*> = k 


x3n,k 


^ K 2 = m = q + — kB 2 



(7.45) 


where intensity and scale for region 2 are to be used. Under the assumption of isotropic turbulence, we 
take the u- and the v-intensities to be the same. 


The equations of Sections 6 and 7 have been programmed in Fortran in code CupBB. Cases that 
exhibit the physical effects of interest typically run 1 to 2 hours on a 200 MHz personal computer. 
However, if high frequencies (significantly beyond the spectrum peak) are calculated, the running time 
can be much longer. 
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SECTION 8 

VERIFICATION AND APPLICATION OF THEORY 


This section provides verification of the new broadband coupling theory and then presents a series of 
calculated power spectra to demonstrate the capability of code CupBB and to demonstrate the 
importance of unsteady coupling and flow turning. Conclusions will be summarized in Section 9. 

Verification 

As described in the background section (Section 2), the new theory is a combination and extension 
of preceding theoretical treatments of tone and broadband noise generation. Reference 1 addressed 
unsteady coupling for tone noise in a 2D environment, with turning included via actuator disks. 
References 2 and 3 addressed broadband noise generation by isolated (uncoupled) cascades. Geometry 
and mean flow were uniform in the spanwise direction but use of a true 3D turbulence spectrum was 
enabled by Glegg’s Wiener-Hopf cascade analysis. In the combined theory of this report, the 
turbulence spectrum, cascade response, and actuator disk response are 3D with respect to unsteady 
perturbations. 

Because broadband coupling has not been treated previously in the literature, there are no test 
cases for verification. However, 3 tests were executed to provide some confidence in the code. First, 
note that the coupling system described in Section 6 must reduce to the tone system of Ref. 1 in the 
case of zero spanwise wavenumber. This has been carefully checked so as to verify (in the special case 
of jt/ = 0, i.e., 2D waves) the coupling, the treatment of actuator disks, and frequency and mode 
scattering. 

The second verification was to run the new code with effects of coupling and flow turning “turned 
off’ and compare results with the older code (BBCascade), which does not treat coupling and flow 
turning. The comparison shown in Figure 16. Ideally, the curves would be identical but there are 
significant differences that must due to differences in treatment of duct boundaries. As can be seen in 
the figure, the differences disappear at high frequencies but amount to about 3 dB at lower frequencies. 
The low frequency increase with the end walls may be a legitimate noise effect associated with noise 
generation in a confined region. Also, it could be that the adaptation of the continuous Liepmann 
turbulence spectrum for the confined environment (as described in conjunction with Figure 15) needs 
to be refined. Finally, note that BBCascade considers a cascade without endwalls and without the 
periodic boundary condition associated with wrapping the rectilinear geometry onto an annulus. 
Without these boundaries, the radial and circumferential wavenumber spectra are continuous. In 
contrast, CupBB enforces the endwall and periodic boundary conditions with the result that discrete 
modes replace the continuous wavenumbers. Modal summation results would approach continuous 
integration at high frequencies. In any case, it is believed that Figure 16 verifies correct behavior of 
the isolated cascade application of CupBB. 

The third verification was to run the new code with all features enabled and check for qualitative 
agreement with test data. Results shown in Figure 17 are again for the same ADP scaled model case 
presented in Section 2, Figure 6. The mode order/frequency plot in Figure 9 still applies. Recall for 
that figure that the turbulence intensity and scale were adjusted to provide a “best fit” since we are 
attempting to represent the fan annulus by a single representative radius via the rectilinear theory. 
Because the modeling with the new theory includes more features, it was necessary to re-adjust the 
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Frequency - Hertz 


Figure 16. Comparison of old and new codes calculation for isolated cascade. BBCascade results 
repeated from Figure 6. Turbulence RMS level was 2.0% of the mean velocity and the integral scale was 
3.5% of the fan mean radius. 

turbulence properties. Intensity was reduced from 2.0% to 1.4% and the integral scale was reduced 
from 3.5% to 2.7% of radius. It was decided to match the downstream data as closely as possible and 
to let the upstream calculation follow, since other sources could be contributing to inlet noise. The 
match to the downstream spectrum shape is excellent, considering that the only free parameter 
governing shape is turbulence scale. This shows the merit of using a 3D turbulence spectrum and a 
3D cascade response function. Earlier schemes based on 2D cascade theory and correlation length 
concepts don’t do this well. 

The significant underprediction in inlet noise should not be considered a failing of this theory. 
Rather it is an indication that turbulent inflow at the stator is an unlikely source of inlet noise, now that 
we have a reasonable model of rotor transmission loss. Recall that, when the isolated cascade theory 
was first brought on line, we were pleased with the upstream/downstream split (in comparison with 
data as shown in Figure 6). This seems naive now that the modeling is more complete. The earlier 
results ignored the effect of the rotor with the result that conclusions regarding upstream noise were 
probably incorrect. For this case, trailing edge noise from the rotor may dominate in the upstream 
direction. 

Regarding the independent mode subsets explained in previous sections, calculations in Figure 17 
were made with 2 sets of 15 frequencies. The first set was half way between the BPF harmonics and 
the second set was interlaced with centers at the BPF harmonics. All of the cut on radial modes were 
included. For circumferential orders, all cut on modes were included plus 2 co-rotating and 2 counter- 
rotating cut off modes for each frequency. (Some cut off modes were included between rotor and 
stator because they may not decay completely before interacting with the adjacent blade row.) Since 
this calculation was very time consuming, short cuts were explored as described in the next section. 
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Frequency - Hertz 


Figure 17. Calculation with fully coupled theory compared with scaled model ADP data from Figure 6. 
Relative to the isolated cascade calculation, the turbulence intensity was reduced from 2.0% to 1 .4% 
and the integral scale was reduced from 3.5% to 2.7% 


Variations in Computational Parameters 

As mentioned above, Figure 17 was computed with NH, the number of frequencies in a coupled 
subset, equal to 15. When NH was reduced by roughly one half (to 8), the computation time dropped 
by a factor of 16. Figure 18 was generated with NH = 8, 12, and 15 to check the effect of reducing 
frequency range. Even though the rotor scatters up and down in frequency, the figure shows that 
calculations are accurate up to whatever frequency is included in the coupled sub-sets. Apparently, 
up-scattering is more important than down-scattering. The remaining calculations is this report are for 
NH = 8. 

An experiment was performed to determine if modal averaging could be used to reduce 
computational effort. Figure 19 was generated by stepping the q index in jumps of qStep and then 
multiplying the results by qStep. It can be seen that skipping every other circumferential mode 
(qStep = 2) has negligible effect at all frequencies. Even skipping in steps of 5 produces good results 
from the peak of the spectrum and up in frequency. Jumping in steps of 10 produces unreliable results. 
This is important information since computation time goes as 1 /(qStep). Despite this benefit, the 
remaining calculations were made with all q’’ s. It also seems likely that radial order could be skipped 
for further computational efficiency but this has not yet been verified. Simply skipping every other 
radial mode and then multiplying by 2 did not provide good results. 
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Figure 18. Effect of varying NH on computed results at lower frequencies. For NH = 8, the curves stop at 
3037 Hz, for NH = 12 at 4555 Hz, and for NH = 15 at 5695. 



100 1000 10000 


Frequency • Hertz 

Figure 19. Effect of modal averaging with qStep = 1,2,5, and 10. For qStep=1 , all circumferential 
modes are included. For qStep=2, every other mode is skipped and the results are doubled, etc. 
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Importance of Coupling and Flow Turning 

Figure 20 was prepared for a direct comparison of the earlier isolated cascade technology and the new 
methodology with all features included. The solid curves, are the data-matched results from Figure 17 
and include all of the new capability in CupBB. The dotted curves are the result of an isolated stator 
calculation with CupBB. It can be seen that the full methodology increases downstream noise by 
about 3 dB at all frequencies and increases the upstream/downstream split (or differential) from about 
5 dB to about 10 dB. 



100 1000 10000 

Frequency - Hertz 


Figure 20. Calculations with CupBB for stator inflow turbulence. Isolated stator calculation compared 
to calculaton with all features included (rotor blades and actuator disks). 


A slight variation on this result is shown in Figure 21. Here we represent the stator with turbulent 
inflow by the cascade plus the stator actuator disk (dotted curves). The solid curves result from adding 
the rotor cascade (with its actuator disk) without changing the mean flow or turbulence at the stator. 
Again the aft noise increased and the upstream/downstream split is increased. 

CupBB treats rotor inflow turbulence as well. The sample calculation in Figure 22 compares the 
previous stator results with noise computed for the rotor with the same normalized intensity and scale 
for the turbulence. Both inlet and aft levels increase about 4 dB indicating greater radiation efficiency 
by the rotor. Note the spectrum shapes for rotor and stator noise are almost the same. Of course, we 
expect average turbulence intensity at the rotor face to be considerably less than at the stator face for 
typical turbofans. Thus, Figure 22 does not imply that inlet noise is dominated by rotor inflow 
turbulence. 
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Figure 21 . Effect of rotor on noise from a stator with turbulent inflow. 
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Figure 22. Comparison of calculations for rotor inflow turbulence with stator inflow turbulence. 
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Finally, in a form parallel to Figure 21, we evaluate the effect of the stator on rotor noise. The dotted 
curves in Figure 23 are for a rotor, represented by a cascade plus actuator disk, with turbulent inflow. 
The solid curves result from adding the stator (with actuator disk) on the downstream side without 
changing the rotor mean flow or inflow turbulence. The effect is very strong and again indicates that 
isolated blade row calculations are too simplistic for reliable fan broadband noise prediction. 



Figure 23. Noise due to rotor inflow turbulence, with and without the effect of the stator. 
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SECTION 9 

CONCLUDING REMARKS 


Earlier methods for predicting turbulent inflow noise from fans treated rotors and stators as 
isolated blade rows in a uniform background flow. This report extends that capability by adding the 
effects of scattering by the adjacent blade row and turning of the mean flow by the rotor and stator. 
Critical elements of the theory include: unsteady acoustic and vortical responses of rectilinear cascades 
via S. A. L. Glegg’s Wiener-Hopf theory, unsteady actuator disks to represent the effect of changes in 
mean flow on perturbation waves, full unsteady vortical and acoustic coupling of the rotor and stator, 
and turbulence represented by its 3D wavenumber spectrum. The modeling requires that geometry and 
mean flow be constant in the spanwise direction. Although this is a limitation, the benefit is that the 
cascade unsteady loading and acoustic/vortical response are mutually consistent and fully 
3 dimensional. 

In developing the coupling analysis, an important principle was discovered: Many modes and 
frequencies are involved in fan broadband noise generation. In the process of reflecting back and forth 
between the rotor and stator, each mode/frequency couples to many other mode/frequencies. However, 
the coupling is limited to independent mode subsets each of which only communicates with its own 
modes and not with the other mode subsets. This principle of independent mode subsets makes the 
coupling problem tractable. If all modes could communicate with each other, computation time would 
be excessive. But, by taking advantage of this principle, all of the calculations presented herein could 
be performed on a personal computer. 

The computer code CupBB, which embodies the broadband theory, can compute noise due to 
turbulent inflow at the rotor or stator or both. Comparison of isolated rotor or stator calculations with 
fully coupled calculations indicate that coupling increases the aft noise by 3 to 4 dB and increases the 
upstream/downstream sound power split from about 5 dB to about 10 dB. 

Based on the limited calculations presented herein, the 2 major conclusions are as follows. 

1. It is unlikely that stator inflow turbulence contributes significantly to inlet noise. 

2. Unsteady coupling augments noise generation to the point that it must be included in any 
attempt at absolute level predictions. 

The analysis presented in this report is based on rectilinear, flat plate theory. Although this has its 
limitations, the 2 conclusions above should be reliable because the analysis is a fully consistent, 3D 
treatment of the unsteady flow. Trends regarding effects of adjacent blade rows and flow turning are 
credible. Thus, in the short term, the analysis can be used to upgrade existing isolated blade row 
prediction methods. In the longer term, the principles discovered here and the coupling method 
derived can be used as guidance in development of unsteady CFD methods that can represent the blade 
response and turning effects more accurately. 

Because of time constraints, conclusions have not yet been drawn regarding the importance of 
mode trapping in fan broadband noise. The physics is modeled in the code but the computer runs 
required for the analysis must be performed in future work. 
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APPENDIX A 

GLEGG’S CASCADE THEORY FOR PRESSURE WAVES ADAPTED TO 
DUCT COORDINATES AND THE n, k INDEX NOTATION 


Section 6 of this report requires forms for the acoustic and vortical responses of stators and rotors to 
excitation by input acoustical and vortical waves. This appendix provides the acoustic response and 
Appendix B provides the vortical response. 

Glegg developed his theory in “cascade coordinates” in which the x axis is aligned with the flow 
direction. In this appendix we adapt his scheme to “duct coordinates,” which have been chosen for the 
coupled cascade application. In duct coordinates the .v axis is aligned with the duct axis, the y axis 
with the tangential direction (positive in the direction of rotor rotation), and the r axis coincides with 
the leading edge of the reference vane (or blade). Also, we replace the usual interblade phase angle, 
common in unsteady cascade analysis, with the n,k index notation commonly used in fan acoustic 
interaction theory. Glegg’s theory accounts for sweep, but, in this report, we deal strictly with unswept 
blade rows. 

First, we will review Glegg’s theory (in cascade coordinates) and then transform his equations to duct 
coordinates. As shown in Figure A-l, geometry is constant in the z direction and the background 
flow is uniform: U=(IT,0,0), as expressed in cascade coordinates. Airfoils are unloaded flat plates. 
Cascade gap, chord, and stagger angle are g, c, and 6 . 
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Figure A-1 . Geometry and flow for 
Glegg's cascade acoustic analysis* 41 . 


The unsteady flow is harmonic in space and time with upwash given by Equation A-l 

M’(x,r) = w 0 e i(Yo*c+ayc+v---o>t) 


(A-l) 


This represents a plane wave that is harmonic in time with frequency (o/ln and upwash complex 
amplitude w a . It is also harmonic in space with x c , y c , and r wavenumbers equal to y a , a, and v. By 
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use of Wiener-Hopf analysis, Glegg derived an equation equivalent to the following for the velocity 
potential of acoustic waves scattered by the cascade in response to the input wave of Equation A-l 


^( X>r) = ± g^£_ £ 0 i[-4(x c -y c d/h)+(<?-2!Tk)y c /h+v:] -icot 




where the upper/lower sign applies to upstream/downstream going waves and 


M =W/a 


P=J 1^ 




= ^ 2 +£ 2 /; 2 tan x e =d/ph A* = kM +rjf 




k 2 =k 2 -{v/pf 


K = co/(ap^) 


Ik ~ ~fk sm Xe ± C( > s \l*e-fk 2 fk =(°- 2nk + KMd )/ s c 


(A-2) 


(A-3) 

(A-4) 

(A-5) 

(A-6) 


and G = y 0 d+ah is the interblade phase angle. D is the Fourier transform of the discontinuity in 

potential across the blade and wakes (in the form of an infinite product) and is the major result of 
Glegg’s derivation. We adhere to Glegg’s notation closely although we use k for the scattering index 
where Glegg used m, and 6 for stagger, where Glegg used %, and a> for radian frequency where Glegg 
used of. Also, the D function of this report is non-dimensional; to obtain the D function of Glegg’s 
report, multiply the non-dimensional version by w 0 c 2 . 

The velocity, pressure, and density perturbations associated with the acoustic wave can be 
obtained from Equation A-2 via 

u = V0 p--p r D<t>/Dt p' = p/a 2 (A-7) 


where p r and a r are the ambient density and speed of sound in a region of the flow denoted by the 
subscript r. This leads to the following form for pressure in the scattered waves 


p ± (x,r) = ± 


inprWgC 2 

ps e 


*=— yl K e~fk 


A* ( x c - y r d / h)+(o-2nk)y c /h+v:] e ~icot ^ ^ 


The formulation above gives the acoustic waves scattered by a cascade for a single planar wave 
input per Equation A-l. Scattering index k runs over an infinite range but, as usual in this kind of 
formulation, only a finite number of waves are cut on (propagate undiminished); the remaining waves 
decay exponentially and, thus, carry no acoustic energy. Cuton is governed by the argument of the 

square root ^K 2 - f 2 ; when the frequency is high enough, the argument is positive and the waves are 

cut on. 


To apply Glegg’s theory in duct coordinates (as used in the main body of this report), we 
transform the exponential in Eq. A-l to 
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e >(k x x+k y y+v:-<or) (A _ 9) 

where x is the axial duct coordinate and the tangential duct coordinate is y=R <p. Furthermore, to deal 
directly with circumferential modes, we recognize that k y y=m(j). This, with (j> =y/R , gives 

k } , = m/R (A- 10) 

(A- 1 1) 

in terms of k x , m, and co . 


and the working form for the exponential in this appendix: 

/ ( k x x+m(p +vz-(dt ) 

+ 4 - I 9 2 

Later on, we will find expressions for s e , £jr , , and \JK~ - fj. 


Since the analysis of this report tracks scattering by the rotor and stator, we need a notation that 
includes the rotor and stator scattering indices. Consistent with the presentation in Section 4, 
circumferential mode order is indexed on q, n, and k via 

m - q + nB\ — kBj (A- 1 2) 

where n is the rotor scattering index, k is the stator scattering index, and q is the mode offset index. 
Axial wavenumber, with our n, k subscripts was given in Eq. 5-10 

M x (~cb + mM y ) + + mM y ) 2 - (m 2 + v 2 ) (A- 13) 



where we have normalized frequency by radius and speed of sound 


and spanwise wavenumber by radius 


„ caR 
(0 = — 

(A- 14) 

a r 


<> 

II 

< 

(A- 15) 


The reader can verify that Eq. A- 13 can also be obtained from Glegg’s wavenumbers by coordinate 
transformation as in the following section. 


Stator Scattering 



The connection between Glegg’s cascade coordinates and our 
duct coordinates is a simple rotation about their common r axis 
as shown at the left. The transformation is given by 
x = x ( . cos 0 2 - y c sin 
y = .v ( . sin $2 + >’ c cos<9 2 

and (A- 1 6) 

x c = x cos 0 2 + y sin 0 2 
y c = -xsin0 2 + >’Cos<9 2 
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These equations are shown for the coordinate transformations but apply to transformation of the 
wavenumbers as well. Note in Equation A-2 that the axial wavenumbers are —X k . Hence, from 
Equation A- 16, 


tfk = ~ k xnk cos 0 2 - j sin d 2 (A- 1 7) 

where we have included an n subscript X^ k to help track scattered waves. This leads to an alternative 
form for the effect of the convective derivative in Eq. A-7 co +WX^ k (times -/) 

co+WX k = o)-W(k^ nk cos 0 2 + — s'm9 2 ) = — — (S + M x k± nk R + M y m) (A- 18) 

R R 

or 

c o + WX k = — t-A * k (A-19) 

A 

where A± k = -cb + M v k± nk R + M v m is the same as defined in Section 5. 

We can also apply the transformation of Eq. A- 16 in the reverse direction to verify the form of the 
tangential wavenumber in the duct system derived from Glegg’s forms. Note in Eq. A-2 that the y c 
wavenumber in the cascade system is X„ k tan 9 2 +(<j - 2nk)/h . Thus, 

os0 2 

or, in terms of vane gap, g 2 

m _ _ <7 — 2 nk 

~R = y ~ 7i 

since h = g cos 9 . This is the well known form from Smith’s theory (Ref. 7). Another expression 
required below is easily derived 

s e =fag 2 (A-22) 

where 

P x =yjl-M 2 x (A-23) 

To deal with the square root in the denominator of Equation A-2, we substitute definitions from 
Eqs. A-3 to A-6 and re-arrange to find 

yR - Tk = + mM ) 2 - ft; (m 2 + v 2 ) (A-24) 

K PPx 

The square root on the right side is the discriminator for cutoff. It appears frequently so we give it a 
special symbol 


77 = - Aj sin 9 2 +( rf k tan 9 2 + 


G -Ink 


(A-20) 


(A-21) 
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Then, 


R nk = ^(- 0 ) + mM y ) 2 -(i 2 (m 2 +v 2 ) 


(A-25) 



(A-26) 


A convenient form for £* can be developed by similar manipulations. It must be treated with 

some care because the definition implies a square root of a square root, where the branch must be 
defined in the case of cutoff. From the definitions in Eqs. A-3 to. A-6, we force the following form 


f» = +/(-(')* ) 2 (A-27) 


Insert the definition of rjf to get 


C* = P -/*)+/*-[ /* 2 s in 2 T 2/, 


sinj e cos^ + 




cos 2 


(A-28) 


and collect terms to find 


Ca = P y/(^e - Tk ) : sin ^ 2 Xe ± 2 fit ~ fk sin Xe fk cos Xe + fk co ^ Xe ( A ' 29 ) 

or simply 

C? = P ( yj^e-fk sin ^ ± /a- cos ^ j ( A - 3 °) 

Again we add the n subscript and non-dimensional ize on R 

(A-3D 


Then, by substituting definitions from above, we find the non-dimensional form to be 




V 


(A.-32) 


In obtaining a working form for the scattered acoustic pressure, we need to distinguish clearly between 
input waves and output waves using our “ n, k ” notation. We will use unprimed indices to indicate 
input waves and primed indices for output waves. The general form for mode order of waves input to 
the stator is 

m = q + nB\ — kBj (A-3 3) 
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where q is the mode offset index, n relates to scattering at the rotor and k relates to scattering at the 
stator. Interblade phase angle of the input waves is related to mode order by 


Thus, in the “ n, k ” notation. 


a =- 


2nm 

~b7 


CJ=^ L (q + nB l -kB 2 ) 

B 2 


(A-34) 


(A-35) 


We temporarily denote the scattering index (of the output waves) as k". In Glegg’s (and Smith's) 
scattering theory, the scattering index appears only in combination with the interblade phase angle (of 
the input wave) as o — Ink’ . We can relate this to the circumferential mode order of the output wave 
via 

O -rr _ ff 

G - Ink * = — [q + nB { - (k + k')B 2 ] = — — (A-36) 

B 2 B 2 


where we have defined 


m” = q + nB l — ( k + k r )B 2 


(A-37) 


We modify our notation for Glegg’s D function to D± kk ' which indicates input mode n,k and output 
mode n,k'. When we put all of this together, the adapted form for the pressure waves becomes 


Av)= I 

^2 R k r =—e 


A 


ruk+k 


f>± 

hn,k+k A 


D f ] f k,k+k* 


k fnJt + * rX+mr 4 ,+v = - <ot 1 


/ 1 


(A-38) 


■n,k+k‘ 


Since the sum runs over all values of k", it can be re-indexed via k' - k + k" with the result 


2 oo . ± y± 4. , 

±/ v Pr a r c 2 ... V ^nk’ h nk' 'nkk' i[k xnk 'X+m<p+v:~cot ] 

P (X,G- 2^ e 

8 2 K k '=- °° K nk' 


(A-39) 


For the duct simulation in the main body of the report, the endwall are at : = 0 and z - h and the r 
dependence of the upwash is given by cosines as in Eqs. 3.10 and 5.7. Of course, cosines are linear 
combinations of Glegg’s upwash form in Eq. A-l 

vv’(x, t) = w Q e Kr ° x ' +ay c- m) cos(i^) (A-40) 

h 


where ji is the radial mode index running from 0 to <». Now Eq. A-39 can simply be modified to read 


p i ^ _ -- rr-rj w „ k Y D nkk' J[k* nk 'x+m<t>-cot] ^, t Pnz 


■±'~ Pr a r c 2 

82 R k ’=- 


R>ik' 


cos( - ) (A-41) 

h 
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This can be written as 


p ± (\,t)= X Pnkk' 

k — oo 


(A-42) 


where p n kk' is the form used in Section 6 to develop the system of scattering equations. 


= +lK P ^ . a . r S 2 Wnk e ^nk‘x+m<p cot] cos(J M^£ ) (A-43) 


Pnkk' 


8 2 


R 


nk' 


We have defined c 2 = c 2 /R and g 2 = g 2/ c 2 • 


Rotor Scattering 

Coupling between rotor and stator is done in the stator reference frame at mode order m = q + 
nB]-kB 2 and frequency C0 o +nBi£2. Thus, the input order and frequency for the rotor are the same as 
for the stator problem. However, to apply the cascade theory, we must transform the input to the rotor 
frame where the mode order stays the same but frequency shifts because of rotation. Mode order m 
results in interblade phase angle 


2nm 

o = 

Consider the kinematic phase of the input wave 

y/ = mcf) - cot 

Transformation to the rotor frame via (p = <p r + Qt gives 

y/ = m<f> r -(co- mQ)t 


(A-44) 


(A-45) 


(A-46) 


The rotor scatters mode order on the n index, so we denote the output wave modes as 

m" = q + (n + n")B l - kB 2 (A-47) 

However, in the rotor- fixed frame, scattering does not change frequency; hence kinematic phase of the 
scattered waves becomes 


yr =m <p r - (co- mQ)t 


Thus, frequency in the rotor frame becomes 


co,. = co- mQ. 


(A-48) 

(A-49) 
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or, in normalized form 


(A-50) 


cb r = cb — mM T 


Axial wavenumber for the rotor can be obtained from Eq. A- 13 by replacing cb by cb, and M y by 


f.± _ 1 

k xnk “ — 


m 1 


M 


x [-d) r + m(M y -M t )] + j[-<b r + m(M y -M T )f- ft 2 (m 2 + v 2 ) 


(A-51) 


However, it can be seen from Eq. A-48 that 


-cb r + m(M y - M t ) = -cb + mM y 


(A-52) 


and axial wavenumber reverts to the form used for the stator in Eq. A- 13 

kxnk = — ~~2 M , + y ) + + +v 2 ) 

RB..L 


(A-53) 


This is as it must be since the wavenumbers are related to physical lengths of the waves that would be 
the same in any coordinate system. A similar approach shows that yjKj - /„ 2 has the same form for 

the stator, i.e., yj Kj - = R nj ,/Rf$fi x where, as before 


R nk = yj(-cb + mM y ) 2 - fi x (m 2 +v 2 ) 


(A-54) 


The derivation for proceeds the same as for the stator down to Eq. A-30: 


^i=P^ K e-fn s in Xe - fn C °S Xe j 


(A-55) 


except that now, of course, Xe applies to rotor geometry and flow. We add the k subscript for clarity 
and normalize on R as follows 




(A-56) 


Substitution of the appropriate definitions leads to the following as the normalized form for the rotor 
Cnk = — J [ R nk sin ± ( + &r ( M y -M T ) ) COS 6 ] J (A-57) 


where again the upper/lower sign goes with the up/downstream-going acoustic wave. Here, 
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tf=l -M 2 x -{M y -M T ) 2 


(A-58) 


is related to the relative Mach number at the rotor. 


As stated above, interblade phase angle for input waves is 

2nm 


a = ■ 


B, 


(A-59) 


where m for the input wave is still m—q+nBj-kB 2 . We temporarily use n" as the rotor scattering 
index. Then the factor in Glegg’s equations containing the scattering index becomes 


O + Inn = — [q + (n + n)B { - kB 2 ] = lKm 


B 


B 


Here we have defined 


m" = q + (n + n")Bi - kB 2 


The pressure in the rotor frame can now be written 


+ , v _ +l£ £r Ml ,, 

P (X, *) — _ M nk 

g\R 


y ^-n+n’,k Yi+n' ,k ^n,k ,n+n' J^ k x.n+n',k x+m +V 
n'=-oe Rfi+n',k 


We re-index this via n'=n+n" and return to the stator frame via 0, = <f> -Qt to find 


p~(x,t) = 


Yin p r a r cf Y A «'A' Cn'k Dnkn A k t n ’k x+m '^ +v: 
w nk 2*i e 


g\R 


n -—oo 


R 


n'k 


(A-60) 

(A-61) 

(A-62) 


(A-63) 


As with the stator scattering equations, A „>/. gives the effect of the convective derivative. This is the 

same in the rotor and stator frames since, by definition, it is the derivative “following a particle.” The 
equivalence can be seen in our notation as follows. The convective derivative in the rotor frame 
produces 

-co r +Uki' k +<y-QR)!j] 


]2j_ 

D,t 


= i 


(A-64) 


and in the stator frame 


D 

~Dt 


= i 


T TJ + T J 

-0) + Uk~ n ' k +V — 


(A-65) 


These are the same since co r =(0- m'Q . . In a form parallel to Eq. A-40 for the stator waves, we write 

oo 

p ± (x,t)= X Pnkn ( A ‘ 66 ) 

n'=—oo 

where 
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(A-67) 



For the upwash in a ducted system we again use Eq. A-40 so that Eq. A-67 becomes 


* i A 

± . _ p r a r Cj ^n’k C D n f. n ' i[k±„'i.x+m$-(ot\ - , 

Pnkri ~ = M nk e COS{——) 

g 1 R nk h 


(A-68) 


This gives the scattering by the rotor represented in the stator reference frame. It is the desired form 
for setting up the scattering equations in Section 6. 
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APPENDIX B 


GLEGG’S CASCADE THEORY FOR VORTICAL WAVES ADAPTED TO 
DUCT COORDINATES AND THE n, k INDEX NOTATION 


This appendix develops equations needed in Section 6 for vortical waves shed by a cascade due to 
excitation by both acoustic and vortical waves. Results are provided in the contexts of both stator and 
rotor. 

Glegg has extended the Wiener-Hopf theory of Ref. 4 to include vortical response (via a personal 
communication). The output velocity perturbation can be written in a series over the scattering index k 
as follows 


V(x,0= X 


(B-l) 


k =— oo 


where x c and y c are Glegg’s cascade coordinates (.v and y in his notation) and y c and a are the 
associated wavenumbers (corresponding to convection). 


*•% 


a = (a - Ink - y c d)/h 


(B-2) 

(B-3) 


and v is the r wavenumber. Glegg’s result for the velocity vector amplitudes is 


V* = ; 2 ^ (y 2 f)VV f, (■ oey c i c +tfic +«vk c ) 


(B-4) 


/t(7c +cr +v*) 


where 


Ce =-(Yc+ v 2 ) 


(B-5) 


When V* is represented in duct coordinates via Eqs. A-16, the resulting x, y, z components are 


_ -niK(y )u 0 c cqs ^ _ ^2 sin e ) 

h(yj +a 2 +v 2 ) 


(B-6) 


v k = 2K,K( yc )w o c (CC7c sin 0 + ^ 2 cos e) 
h(y c +a +v ) 


(B-7) 


-2niK(y c )w c 
w k = 771 — ^-V(«v) 
h(y^ +a +v ) 


(B-8) 
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Stator Waves 

For the vortical waves, we apply the strategy already developed in Appendix A for pressure waves. 
Mode order of the input waves is 


m = q + «S] — kBj 


with associated interblade phase angle 


<7 =- 


Inm 

~b7 


(B-9) 

(B-10) 


Since the stator scatters on the k index, \fre write the mode order of the scattered waves as 


m” = q + — ( k + k”)B 2 


(B-l 1) 


Then the y c wavenumber can be expressed in terms of duct coordinates wavenumbers via Eq. A- 16 
with the result 

(B-l 2) 


a, 


n,k+k" = "*v3 »,*+*' sin$ 2 + ^COS0 2 


Now the velocity perturbations controlling upwash (in the duct coordinate system) are 
-2niK(y c ) 


and 

v(*,0 = 


2 

V O-nJ.+k’Yc co *®2 — Cc s i n $2 i(k^„ k+k '.x+m'^+v:- 0 )t) 
, / W "k 2- 2 2 2 C 

h i c r=-~ Yt+Kk+k'+v 


-2KiK(y e ) f. 
h/c "* ^ 


&n,k+k*Yc s * n ^2 cos ^2 Kk x $ n £ + k*x+rn 1 '<l>+vz--C)t) 
o Yc+ a ik +k '+ v2 


(B-l 3) 


(B-14) 


Eqs. B-l 3 and B-14 are in a mixed notation, using wavenumbers from both the duct system and the 
cascade system. This turns out to be convenient for coding. 

Again, we shift the k" index via k' = k + k" to arrive at the final forms for vortical waves shed by the 
stator 


and 


where 


«(X, 0 = X u nkk’ 

k'—- OO 


’(X,0 = X V nkk’ 
k 00 

,1 . 


(B-15) 


(B-l 6) 


_ -2niK(y c ) cos0 2 Cr sin0 2 J(*v3n*'*+"'V+v , --®0 ^ ^ 

u nkk r; 2 2 2 M «A- c 

h ! C Yc +a nk' +V 


_ -2mK(y c ) a nk >y c sm0 2 + £ r cos6) 2 ... j(k r , nl ,x+ m ’<t>+v-.-Q>o 

v nkk ~ 77 2 2 2 

h/c Yc +a nk’ +V 


(B-l 8) 
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As the notation implies, these represent scattering by the stator of an n,k input wave into an n,k’ 
output wave. 

Rotor Waves 

Here we find forms analogous to Eqs. B-17 and B-18 for scattering by the rotor. Input waves have 
mode order m = q + nBi - kB 2 and frequency co in the stator coordinate system. Mode order m 
results in interblade phase angle 


<7 =- 


2nm 

~bT 


(B-19) 


Consider the kinematic phase of the input wave 

if/ = nup - cot 


(B-20) 


Transformation to the rotor frame via (j) = (j), ■ + Qt gives 

y/ — m(j> r - (co- mQ.)t (B-2 1 ) 

The rotor scatters on the n index, so we denote the output wave mode orders as 

m'=q + {n + n”)B y -kB 2 (B-22) 

However, in the rotor-fixed frame, the frequency does not change so that the kinematic phase of the 
scattered waves becomes 


if/" = — (co — mQ)t 


Then the a wavenumber is 


a 


n+n’,k ^x3n+n',k Sin0| + ^ COS0] 


(B-23) 


(B-24) 


and the relevant velocity perturbations (in the rotor frame) are 


u(x t) = (Yc ) w ^ ^ a n+rt',k Is. C0S ^1 C S ^ n ^l J(kxin+n,k -V+V :+!//') 


and 


v(x,0 = 


V ' ’ nK AmJ 2 2 2 

c Yc + <W,*+ V 


-2 niK(y c ) ^ a n+n’,k Yc sin 0 l + C? COS0 1 i(k x3n+n ' k x+v: + V ’) 

w nk 2^ : ~ " e 


hjc 


n =-oo 


Yc + a i+n\k +V 2 


(B-26) 


Transformation back to the stator frame yields 

y/' = m'(p - (G) + m'£l - m£l)t 


(B-27) 
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We shift the n" index via n' = n + n" to arrive at the final forms for vortical waves shed by the rotor. 
Frequency becomes 6) r = (O + {ri - n)Q and 


and 


oo 

w(x,0= X *W 

t]'=-oo 


oo 

'’(x»0 ~ v nkn' 

n'=—°° 


(B-28) 


(B-29) 


where 
u nkn = 


_ 2^T iK(y c ') Ct 2^ COS^I sin^, i{k x $ n 'i.x+m'<l>+v:— [<w+(n / — 


hjc 


2 2 2 

rt+apk +V^ 


'n/t 


(B-30) 


_ 2,JtiK(fy c ') a„ y r sin0j + ^ COS0, ^d^rin'i-x+mV+v;— [ct)+(fl ,— n)B]£2]f) 

l nkn' ~ , / 2 2 2 w nk e 

h / c K+rik+v 2 


(B-31) 


And here the notation indicates scattering by the rotor of an n, k input wave into an n', k output wave. 
These forms are adapted to the coupling analysis in the main body of the report. 
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APPENDIX C 

ACTUATOR DISK THEORY 


In Section 6, equations were derived for the rotor and stator as scattering elements. For example, the 
stator element consists of 2 sub-elements: the vane cascade in uniform flow and an actuator disk whose 
role is to account for changes in the mean flow that affect the acoustic environment. Section 6 gave 
the derivation for the cascade scattering and showed how the cascade is coupled to the actuator disk; 
this appendix provides the derivation of the actuator disk equations for application in Section 6. 

The actuator disks in the scattering theory account for turning of the mean flow at the rotor and stator 
through conservation of mass and momentum. In the conservation equations, the flow quantities are 
written as sums of steady and perturbation (or unsteady) parts. When the steady and unsteady parts are 
separated, the steady flow is considered to be prescribed or known from a separate aerodynamic 
analysis. The steady flow provides coefficients for the unsteady equations so that jumps in the mean 
flow cause jumps in the perturbation flow. The jumps are effectively modal reflection and 
transmission coefficients that can be found by inverting a linear system. 

There are 4 wave types: 

T=1 : Upstream-going pressure waves 

T=2: Downstream-going pressure waves 

T=3: Vortical (downstream-going) waves 

T=4: Vortical (downstream-going) waves 


and 4 conservation equations: 

Conservation of mass 
Conservation of axial momentum 
Conservation of tangential momentum 
Conservation of radial momentum 

These are satisfied on a mode-by-mode basis using the standard wave set from the main text. Each of 
the 4 wave types as input scatters out 4 wave types. Thus, whereas the cascades scatter on mode order 
and wave type, the actuator disks scatter only on wave type. The actuator disks do not scatter on mode 
order because they have no variation in geometry in the tangential direction. 

In the following sections, we apply the conservation equations, one at a time, by linearizing them and 
applying the standard wave sets from Section 5. These lead to 4 by 4 matrix systems for each 
circumferential mode order (each n, k combination) that are solved by inversion. 


CONSERVATION OF MASS 

For conservation of mass, p u is matched on both sides of the actuator disk: 

(pu) a =(pu) b (C-l) 

We write this in terms of steady and unsteady parts (with the tilde’s): 
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(Pa + Pa a +«*) = (Pb + Pb^b + "t>) 


(C-2) 


and focus on the first order unsteady terms, which become 


. U a . „ U h . 

Pa w a + — Pa ~ Pb u b + ~Pb 


a. 


a b 


7 

after applying p = pa“ . Now, we define 


f'-' _ 

/7 


- - 

Pa 11 a 2 P * 7 

a„ 


(C-3) 


(C-4) 


for matching at the actuator disk with a similar form C b for region b. The normalization via a 0 !p 0 is 
for convenience. a Q and p a are reference values for the sound speed and ambient pressure which will 
drop out of the formulas later. 

We interpret Eq. C-4 to apply on a mode-by-mode basis and express the perturbation quantities as the 
sum over all wave types. For example, the contribution of the upstream-going pressure (Type 1) in 

region a is denoted by the coefficient Aj 7 , which is stands for the defining pressure component 

Aj° (n,k) in Eq. 6.1. The contribution from the associated axial velocity component in Eq. C-4 is 

proportional to Af via the second line of Eq. 5. 1 1. 


For p = 0: 
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(C-5) 


where the first line comes from the upstream-going pressure wave, the second from the downstream- 
going pressure wave, and the third from the vorticity wave. Because we are matching on a mode-by- 
mode basis, we have dropped the n,k subscripts, the exponential, and the cosine, which are common to 
all terms. 
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For fi> 0: 
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(C-6) 


All of the matching is at .v = 0; reference planes are shifted later via exponentials. We compress the 
notation and write Eq. C-4 as 


where 


For jj = 0 


and for p> 0 


C a = C l a A \ + C 2a A 2 + C 3a A 3 
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(C-7) 


(C-8) 


(C-9) 


(C-10) 


where y the ratio of specific heats, entered from yp a = p a aj . The same argument leads to the 
conserved quantity on the b side of the actuator disk (for p > 0) 

Cb = C\b A \ +C 2b A 2 + C 3b A 3 + C 4b A 4 (C-ll) 


Now we equate C a and Cb- Furthermore, we identify the 4 waves approaching the actuator disk from 
each side as input waves and place them on the right hand side of the equation and the 4 waves leaving 
the actuator disk as scattered waves and place them on the left side. The input waves are the upstream- 
going pressure wave on the b side and the downstream-going pressure wave and the vorticity waves on 
the a side. The scattered waves are on the upstream-going wave on the a side and the downstream- 
going pressure and vorticity waves on the b side. 
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(C- 12) 


C\a A \ -CibA 2 ~ C 3 b A 3 ~^ 4 b A 4 = C \b A l ~ C 2 a A 2 ~ C 3 a A 3 “Q« 


This (or the corresponding equation for p = 0) is the first of the actuator disk equations and represents 
conservation of mass. 


CONSERVATION OF AXIAL MOMENTUM 

We treat conservation of axial momentum in similar fashion. The unlinearized conservation equation is 

Lsina m + (p + pu 2 ) a =(p + pu 2 ) b (C-13) 

where the first term is the axial component of loading on the disk. The critical step here is to separate 
the total loading on the blade row into the steady part, which is handled via the actuator disk, and the 
unsteady part, which is handled via unsteady cascade theory. Then, when we express the variables in 
Eq. C-13 in terms of steady and unsteady parts, the unsteady loading (on the actuator disk) is by 
definition equal to zero and the first order perturbation equation is 

(l + (M a x ) 2 )p a +2p a U a u a = (l + (M b x ) 2 )p b + 2p b U b ii b (C-14) 


We define 


= —[(! + ( M x ) 2 )Pa + 2 Pa U a ii a 

Po l 


and expand it in the 4 wave types, as before 
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for p = 0 


and for p > 0 
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(C-15) 


(C-16) 


(C-17) 


(C-18) 


(C-19) 
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The equations for the b side of the actuator disk are the same with b's substituted for the a' s. By using 
the scheme described above for identifying source waves and scattered waves, the equation for 
conservation of axial momentum becomes 

F ia A° - h,A - F,A - Ft. A = FuA - F 2 a4 -FlaA- FtaK (C-20) 


CONSERVATION OF TANGENTIAL MOMENTUM 

The unlinearized equation for transverse (y direction) momentum is 

-Lcosa m + {puv) a =(puv) b (C-21) 

Using the same arguments as in the discussion of axial momentum, the first order conserved quantity 
can be written 


G a = —\pa U a^a + Pa V a l ~ l a + M x M y p a ~^ 
Po L 

or, in terms of the 4 wave types, as 


where 


for {1 = 0 
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(C-22) 


(C-23) 


(C-24) 


(C-25) 


and for {i> 0 


Po a a L J 
Po a a l ' J 


(C-26) 


and the conservation equation for tangential momentum has the same form as the other 2 conservation 
equations 


G la A l ~ G 2b A 2 - G 3b A 3 ~ G 4b ^4 ~ G l b A \ ~ G 2a A 2 ~ G 3a A 3 ~ G 4a A 4 (C-27) 
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CONSERVATION OF RADIAL MOMENTUM 

This, the 4th conservation equation, is only applied for fi > 0, where there are 4 wave types. The 
unlinearized equation for transverse (y direction) momentum is 

(puw) a = (puw) b (C-28) 

Using the same arguments as in the discussion of axial momentum, the first order conserved quantity 
can be written p a U a w a . We define 


or, in terms of the 4 wave types, as 


H a [ Pa G a a J 

Po 


(C-29) 


Ha = H la A? + H 2a A% + H 3a A% + H 4a A% 


(C-30) 


where 






K Kk 




Po a a 

H 4a = 

Po ® a 



(C-31) 


and the conservation equation for radial momentum has the same form as the other 3 conservation 
equations 


H\atf - - «3/,4 - «4f A = - H 2a A f -H 2a A° - H 4a A% (C-32) 


SOLUTION OF LINEAR SYSTEM 

By “solution” here, we mean finding the scattered waves for specified input waves. To this end we 
assemble Eqs. C-12, C-20, C-27, C-32 into matrix form 
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Solution to this equation is 
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(C-34) 


(C-35) 


APPLICATION OF SOLUTION 

The objective of this analysis is to create rotor and stator acoustic elements that can be represented in 
the coupling system of Figure 13. For example, when the stator cascade (with uniform flow 
corresponding to region 2 ) is combined with an actuator disk at the trailing edge to turn the flow from 
the 9 direction to axial, the result is the stator acoustic element. This is done by using Eqs. 6-13 and 
6-14, which amounts to 4 equations for the input/output behavior of the actuator disk, and a similar set 
of 4 equations for the input output behavior of the stator cascade. Setting the output of the cascade 
(on the downstream side) to the input of the actuator disk (on the upstream side) and the output of the 
actuator disk to the input of cascade permits us to eliminate 4 equations. The result is Eq. 6-22 which 
represents the input/output characteristics of the combined stator element. 
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APPENDIX D 
LIST OF SYMBOLS 


a n a o 
Cj ,C 2 

8 

h 

i 

k 

m 

n 

P 

<7 

r 

t 


u, v, w 
X, y, : 


x c , y c , z 
w 


speed of sound in region r , speed of sound on standard day at sea level 
chord of rotor and stator cascades 
cascade gap 

duct height in main part of report, gap perpendicular to blades in Appendices A and B 

index for counting frequency from one mode sub-set to another 

scattering index for stator 

circumferential mode order 

scattering index for rotor 

pressure 

mode offset index 

subscript (or superscript) for region (see Section 5). Also denotes rotor, 
time 

axial, tangential, and radial components of perturbation velocity 
coordinates in “duct system” (x axis aligned with axis of rotor rotation) 
coordinates in “cascade system” ( x c axis aligned with airfoil chord and mean flow) 
upwash velocity component 


A r = a r / a v 

A^(n,k) modal coefficient in state vector. See Section 6. T is wave type, r is region, and 
n and k are the rotor and stator scattering indices. 

Bj (n, k ) modal coefficient for source vector (related to turbulence in Section 7). 

Bi ,Bz Number of blades in rotor, vanes in stator 
D Glegg’s potential jump. See Eq. A-2 of this report and Ref. 4. 

E Cutoff discriminator. See Eq. 7.10. 

F n k Ratio of modal sound power to modal sound pressure, see Eq. 7.9. 

H d h/R , annulus height/effective radius 

I acoustic intensity 

K Glegg’s vorticity factor. See Eq. B-4 and Ref. 4. 

L turbulence scale/)? 

M x , M y axial and tangential Mach numbers of mean flow, U/a r and Vi a r 
Pj() See Section 5 


P, =Pr/Po 

R “effective radius” of fan; used for scaling 

/ __ 

Sj'ij- ( ) scattering coefficient, see Section 5 
T wave type. See Section 5. 

U axial component of mean velocity 

V tangential component of mean velocity 

W mean velocity 


A state vector, see Section 6 

B source vector, see Section 6 

S scattering/coupling matrix, see Section 6 
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a chordwise wavenumber in cascade coordinate system 

/ 3 Vl-M 2 

Px 

7 ratio of specific heats, 1 .4 for air 

8jj Kroneker delta, = 1 for / = j and = 0 for /' * j 

8{ ) Dirac delta or impulse function 

= 1 for jj = 0 and = Vz for ju > 0 
= 2 for /t = 0 and = 1 for ju > 0 

i p = y/R, tangential angle in duct coordinates 

H radial mode order 

v radial wavenumber = /jk /Hp 

X result of convective derivative acting on exponentials, see Eq. 3-7 

a interblade phase angle 

6 stagger angle 

p density 

x one half the time for the block of flow under consideration to pass through the 

cascade, see discussion in conjunction with Eqs. 7-25 to 7-29 
ft) 2 n times frequency 

ft ) 0 offset frequency 

Q angular speed of rotor 

A XR / a r 


superscripts 

( ) r region =1 upstream of rotor, =2 between rotor and stator, =3 downstream of stator 
( ) complex conjugate 

subscripts 

( )t denotes wave type, see Section 5 

overbars and hats 
( ) normalization by chord 

( ) normalization by R 
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